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Abstract 

This document provides a summary of the discussions during the recent joint QCD Town Meeting 
at Temple University of the status of and future plans for the research program of the relativistic 
heavy-ion community. A list of compelling questions is formulated, and a number of recommendations 
outlining the greatest research opportunities and detailing the research priorities of the heavy-ion 
community, voted on and unanimously approved at the Town Meeting, are presented. They are 
supported by a broad discussion of the underlying physics and its relation to other subfields. Areas 
of overlapping interests with the “QCD and Hadron Structure" ("cold QCD”) subcommunity, in 
particular the recommendation for the future construction of an Electron-Ion Collider, are emphasized. 
The agenda of activities of the “hot QCD” subcommunity at the Town Meeting is attached. 
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1 Executive Summary 

On September 13-15, 2014, the QCD Community within the APS Division of Nuclear Physics met 
at Temple University, Philadelphia, for a three-day Town Meeting to discuss the status and future 
priorities of their research program, in preparation for a new Long Range Plan for Nuclear Physics to 
be written and submitted to the NSF and DOE Nuclear Science Advisory Committee in 2015. The 
U.S. Nuclear Physics QCD Community consists of two, partially overlapping, subcommunities whose 
activities focus on “The Phases of QCD Matter” (a.k.a. as “hot QCD" community) and “QCD and 
Hadron Structure" (a.k.a. as "cold QCD” community), respectively. Their joint Town Meeting featured 
a daylong joint session (spread over two days, from 4pm on Saturday to 4 pm on Sunday), with separate 
sessions of the “hot" and “cold” QCD subcommunities scheduled in parallel for the rest of the time. 
The Town Meeting program schedule, including links to the slides of all presentations, can be found 
at https ://indico .bnl. gov/conf erenceDisplay ,py?confId=857. With 244 registrants (108 from 
“hot QCD”, 136 from “cold QCD"), participation by the community in the 2014 joint QCD Town Meeting 
was about 10% stronger than for the corresponding meeting in early 2007 at Rutgers University, reflecting 
the health and strength of the U.S. QCD Community in Nuclear Physics. Strong community interest in 
the planning for the next decade was also reflected in a total of 49 short presentations submitted and 
delivered by participants, in addition to the 47 invited talks (15 and 19 in the “hot” and “cold” QCD 
parallel sessions, respectively, and 13 in the joint session) that were solicited, with detailed charges to the 
speakers to ensure a well-rounded program, by the Town Meeting conveners (Haiyan Gao, Ulrich Heinz, 
Craig Roberts, and Paul Sorensen). In this Executive Summary we list the four recommendations that 
were voted on and approved by the “hot QCD” community during the parallel and joint sessions. 

Members of the “hot QCD” community reported on recent progress in the field including the development 
of dynamical models accurate enough to provide quantitative access to emergent phenomena in Hot QCD. 
Heavy-ion data were shown to prefer an equation of state that agrees with that extracted from lattice 
QCD calculations and it was firmly established, through improvements in the determination of transport 
properties, that the hottest matter ever made by humans is also the most perfect liquid known in nature. 
Following on the significant progress made in heavy-ion studies probing zero baryon chemical potential, 
the Hot QCD community is seeking to 1) measure the temperature and chemical potential dependence of 
transport properties especially near the phase boundary, 2) explore the phase structure of the nuclear 
matter phase diagram, 3) probe the microscopic picture of the perfect liquid, and 4) image the high 
density gluon fields of the incoming nuclei and study their fluctuation spectrum. 

These goals, which can be addressed with a targeted experimental program and continued strong support 
for theory, led to the unanimous approval of the "hot QCD” community’s first and highest priority 
recommendation: 

Recommendation #1: 

The discoveries of the past decade have posed or sharpened questions that are central to 
understanding the nature, structure, and origin of the hottest liquid form of matter that the 
universe has ever seen. As our highest priority we recommend a program to complete the search 
for the critical point in the QCD phase diagram and to exploit the newly realized potential of 
exploring the QGP’s structure at multiple length scales with jets at RHIC and LHC energies. 
This requires 

• implementation of new capabilities of the RHIC facility (a state-of-the-art jet detector such 
as sPHENIX and luminosity upgrades for running at low energies) needed to complete its 
scientific mission, 

• continued strong U.S. participation in the LHC heavy-ion program, and 

• strong investment in a broad range of theoretical efforts employing various analytical and 
computational methods. 
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A high-energy, high-luminosity, polarized Electron Ion Collider is considered a requirement to achieve a 
full understanding within QCD of the structure and properties of hadrons and nuclei. The construction 
of an EIC as a means to usher in a new era of precision measurements was unanimously approved as 
the highest priority for future new construction (after the completion of FRIB) by the entire U.S. QCD 
Community in Nuclear Physics: 

Recommendation #2: 

A high luminosity, high-energy polarized Electron Ion Collider (EIC) is the U.S. QCD commu¬ 
nity’s highest priority for future construction after FRIB. 

Recommendation #3 is an amended version of a recommendation passed by overwhelming majority in 
the joint session to “endorse the new initiatives and investments proposed in the Recommendation and 
Request received from the Computational Nuclear Physics Town Meeting.” The amended version listed 
here was passed unanimously in a later vote taken during the separate voting session of the “hot QCD" 
community: 

Recommendation #3: 

We endorse the new initiatives and investments proposed in the Recommendation and Request 
received from the Computational Nuclear Physics Town Meeting, at a level to be determined 
by the requested NSAC subcommittee. In addition, we recommend new funding to expand the 
successful “Topical Collaborations in Nuclear Theory” program initiated in the last Long Range 
Plan of 2007, to a level of at least one new Topical Collaboration per year. 

The final recommendation regarding education and innovation was passed unanimously in identical form 
by both subcommunities in their parallel voting sessions: 

Recommendation #4: 

The QCD community endorses and supports the conclusions from the Education and Innovation 
Town Meeting. 

The rest of this document provides supporting arguments for these recommendations. It does not simply 
summarize the presentations and discussions at the Town Meeting, but tries to provide a broader context 
and to present, in the words of the authors listed on the cover, a coherent picture of the status and future 
of the field. In formulating this text, the authors made extensive use not only of the presentations at the 
Town Meeting, but also of additional material collected for the community white papers by the "Hot 
QCD” and "Electron Ion Collider” communities 00 whose scope is significantly more comprehensive 
and to which the interested reader is directed for further details and additional projects. Section [2] offers 
an overview of the status, recent achievements, major open questions and future goals of research on the 
“Phases of QCD Matter", and describes how the above recommendations arise from such an overarching 
view. Section [3] discusses in greater depth the science behind these future goals and what is needed to 
successfully address the open questions listed in Section [2] Section [4] provides brief descriptions of the 
science at and the design options for an Electron Ion Collider. Like our Recommendation #2 above, 
the material provided in that section, although slightly rearranged, is identical in physical content with 
a corresponding section in the “cold QCD" summary of the joint QCD Town Meeting, to reflect the 
unanimous and strong support for the EIC project by the entire U.S. QCD Community. References and 
an Appendix with the Town Meeting schedule are found at the end of this document. 
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2 Overview and Recommendations 


2.1 Where we stand: recent insights and open questions 

The bulk of the mass of the visible matter in the universe comes from energy stored in the strong interaction 
between its fundamental constituents. At a basic level, the strong interaction is well understood and 
described by Quantum Chromodynamics (QCD). Strongly interacting multi-particle systems feature 
numerous emergent phenomena that are difficult to predict from the underlying QCD theory, just like in 
condensed matter and atomic systems where the interactions are controlled by QED theory. 

We now know that, as the universe evolved to its present state, strongly interacting matter existed in 
at least two distinct forms, a quark-gluon plasma (QGP) phase containing deconfined color charges 
that filled the universe homogeneously during its first few microseconds, and a clustered phase in which 
colored degrees of freedom are permanently confined into color-neutral objects called hadrons which 
make up the nuclei of today's atoms. Additional forms of strongly interacting matter that include a color 
superconducting phase may still exist in the cores of some compact stars or may have briefly existed in 
others before they collapsed into black holes. QCD matter is expected to have a complex phase diagram, 
possibly including one or more critical points. 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory and the Large Hadron 
Collider (LHC) at CERN enable us to study strongly interacting matter at extreme temperatures in the 
laboratory. Conjectured over 35 years ago, the existence of the QGP phase was established unambiguously 
by experiments at RHIC which also made the surprising discovery that the QGP is strongly coupled 
and behaves like an almost perfect liquid. The successful implementation of luminosity and detector 
upgrades at the RHIC facility, recommended in the last Long Range Plan, and the beginning of the 
higher-energy heavy-ion program at the LHC in 2010, made additional studies possible that have begun 
to yield precise estimates for some of its key properties. Decisive theoretical progress at several fronts has 
been instrumental for these achievements. In response to a recommendation in the last Nuclear Physics 
Long Range Plan, the DOE made a dedicated investment into three Topical Collaborations in Nuclear 
Theory, among them the highly successful JET Collaboration. Both DOE and NSF have provided strong 
support for high performance computing. Facilitated by these initiatives, a standard framework has been 
developed and implemented that describes, with quantitative predictive power, the complete dynamical 
evolution of and the interaction of penetrating diagnostic probes with the expanding QCD matter created 
in heavy-ion collisions. 

The unparalleled flexibility of the RHIC facility to collide atomic nuclei of different sizes over a wide 
range of energies, complemented by p+p, p+Pb and Pb+Pb collisions at the LHC with about 15 times 
the top RHIC energy, provides the experimental leverage necessary to clarify the nature of QCD matter. 
New discoveries made over the past decade have sharpened some questions and posed several new ones 
that address the core of our understanding of the nature, structure and origin of the QGP liquid. These 
questions frame our research program for the coming decade. To address them requires, in the short 
term, a suite of facility and detector upgrades at RHIC and the LHC and a series of new experiments that 
exploit these upgrades. In the long term they necessitate the construction of an Electron Ion Collider 
(EIC). The questions, in part motivated by recent discoveries summarized in the following subsection, are 
listed here and expanded upon in the science sections of this document: 

• What are the transport properties of the QGP? How do they change when the plasma is 
heated or doped with excess quarks? 

• How do the collective properties of the QGP liquid, one of the most strongly coupled 
forms of matter now known, emerge from the interactions among the individual quarks 
and gluons that we know must be visible if the liquid is probed with sufficiently high 
resolution? 
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• What is the precise nature of the initial state from which this liquid forms, and how does 
it reach approximate local thermal equilibrium in the short time and rapidly expanding 
environment provided by heavy-ion collisions? 

• Can dense systems of quarks and gluons act like strongly coupled liquids without ther- 
malizing? Does the Color Glass Condensate state that manifests itself when fast-moving 
atomic nuclei are probed at very small longitudinal momentum fraction exhibit collective 
behavior? 

• What is the structure of the QCD phase diagram? Does it, like that of water, feature a 
critical end point separating a line of first-order phase transitions at large baryon density 
from the rapid but continuous crossover found by lattice QCD at low baryon density? 

• How does the observed structure of the QGP change when it is probed at different length 
scales, with photons, jets, and heavy quark flavors? What is the shortest length scale on 
which the plasma liquid looks liquid-like? 

• What is the smallest size and density of a droplet of QCD matter that behaves like a 
liquid? 


2.2 Accomplishments and future goals 

We briefly summarize some recent achievements that play decisive roles in defining the research priorities 
for the coming decade. Additional accomplishments and related future developments are described in the 
science sections of this document. 


Flow fluctuations and correlations 

One of the most important recent breakthroughs has been the realization that each collision event (“Little 
Bang”) exhibits its own unique flow pattern, with measurable strength up to the 7th or 8th multipole. 
It was recently discovered that these ripples in the near-perfect QGP liquid bring information about 
initial-state nucleonic and, likely, sub-nucleonic gluon fluctuations into the final state. This has opened 
new possibilities to study the dense gluon fields and their quantum fluctuations in the colliding nuclei via 
correlations between final state particles. Precise measurements of the complete set of measurable flow 
coefficients, and their event-by-event fluctuations and correlations in both magnitude and flow angle, for 
a number of identified hadron species will make it possible (in conjunction with the following item) to 
map the transverse and longitudinal spatial dependence of the initial gluon fluctuation spectrum. This 
will provide a test for QCD calculations in a high gluon density regime. 


Quantifying the fluidity of quark-gluon plasma 

Precise measurements of the flow coefficients up to the 5th and 6th multipole for charged hadrons, and 
of elliptic and triangular flow for several identified hadron species, have made it possible to determine the 
shear viscosity to entropy density ratio, 77/s, of the QGP to within less than a factor of 2 . This was only 
possible with recently converged results from lattice QCD for the QCD equation of state and due to the 
development of a comprehensive and sophisticated theoretical model for the dynamical evolution of the 
collision, with relativistic viscous fluid dynamics to describe the QGP liquid at its core. A simultaneous 
analysis of RHIC and LHC data has provided first evidence for an increase of (77/s )qgp with rising 
temperature. The present uncertainty on ( t?/s)qgp (conservatively estimated as ± 50 %) is not so much 
limited by the quality of the available experimental data (which would be sufficient to determine the 
viscosity with relative precision of 5-10%) as by the incompleteness of the presently available set of 


measurements. It has recently been understood (see Sec. 3.1) that measurement of a complete set of 
correlation functions between the momentum-dependent anisotropic flow coefficients and their associated 
flow angles, for several particle species covering a wide range of masses, can not only tightly constrain the 
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spatial dependence of the initial gluon fluctuation spectrum but will also result in a large improvement in 
our knowledge of the QGP transport coefficients. 

The recent unexpected discovery of collective, anisotropic-flow signatures in p+Pb collisions at the LHC 
and in d+Au and 3 He+Au collisions at RHIC suggests that similar signatures seen in very-high-multiplicity 
p+p collisions might also be of collective origin. How collectivity develops in such small systems cries out 
for explanation. The additional running this year of p+Au and p+AI at RHIC augmented the d+Au and 
3 He+Au data and added a crucial comparison at the same collision energy. The inescapable last question 
on the list above can only be answered systematically by exploiting RHIC’s flexibility to collide atomic 
nuclei of any size over a wide range of energies. 

Color opacity and temperature evolution of QCD matter from hard and penetrating probes 

Parton energy loss and jet quenching: Hard probes yield information about how energetic partons 
diffuse in transverse momentum space and lose energy as they slice through strongly coupled QGP. 
State-of-the-art analysis of the energy loss of leading hadrons in jets, together with significant recent 
advances in theoretical modeling, both within perturbative QCD and by introducing insights from strong 
coupling calculations, have increased the precision of our knowledge of the transverse momentum diffusion 
parameter, g/T 3 , by about an order of magnitude, to within a factor 2-3. The parton mass dependence 
of jet modification and energy loss in heavy flavor jets will make it possible to separately quantify the 
contributions from different energy-loss mechanisms. More generally, detailed studies of jet modification 
by strongly coupled plasma over a wide range of angular and energy scales will connect its macroscopic 
hydrodynamic description to a microscopic description in terms of quarks and gluons. As such, these 
jet measurements will provide unique microscopic tools to move closer to a fundamental understanding 
of how a strongly coupled liquid can arise in an asymptotically free gauge theory. Such measurements 
require high luminosity operation and new instrumentation at RHIC and the LHC, and their quantitative 
interpretation rests on further development of theoretical tools for a direct comparison of calculations to 
the data. 

Quarkonium thermometry: Precise and systematic measurements of quarkonium production can 
determine the screening length of the static QCD force in a QGP. Screening effects are expected to 
be easier to discern in bottomonium production, due to the absence of diluting effects from bottom 
quark recombination. An initial observation of a significant suppression of the three T states in Pb+Pb 
collisions at the LHC, indicating a sequential suppression pattern, was recently made by CMS; a precise 
measurement will be available from the LHC by 2023. Low resolution and low statistics measurements 
have been made by PHENIX and STAR and are consistent with the suppression of higher T states but 
the different states could not be individually resolved. With increased beam luminosity and the newly 
installed Muon Telescope Detector the statistical significance of STAR'S measurements and its ability to 
separate the different T states will be improved over the next several runs, and by 2021 the sPHENIX 
experiment will be making T measurements at RHIC with both excellent mass resolution and much better 
statistical precision. The combination of these data sets at quite different initial temperatures, together 
with the different sizes of the three T states which provide measurements at three different length scales, 
will provide strong constraints on the screening length in hot QCD matter. 

Electromagnetic probes: Electromagnetic radiation from the Little Bangs integrates over the electro¬ 
magnetic spectral function of hot QCD matter as it changes with position and time. This provides 
information on the temperature evolution of the expanding fireball and opens a direct window on how 
the degrees of freedom in the vector channel change with temperature. Pioneering measurements by 
PHENIX have recently been augmented by first results from Pb+Pb collisions at the LHC and dilep¬ 
ton measurements at several collision energies by the STAR collaboration. These measurements were 
central in proving that the temperature achieved in heavy-ion collisions is the hottest ever man-made 
temperature. More precise future determinations of the low-mass dilepton spectrum are expected to 
lead to an improved understanding of chiral symmetry restoration at high temperature. Total yield and 
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spectral slope measurements in the RHIC Beam Energy Scan (BES) program will help to quantitatively 
determine the changing fireball lifetime and temperature history at decreasing collision energy. Recently, 
unexpectedly large elliptic flows of direct photons measured by PHENIX and ALICE have presented the 
theory community with a puzzle. Its resolution requires more precise future measurements of the yields, 
slopes and anisotropic flow coefficients of direct photons at RHIC and LHC collisions. 

Mapping the QCD phase diagram 

Theoretical models suggest a phase diagram for QCD matter that rivals that of water in complexity. It is 
the only experimentally accessible phase diagram of matter that is controlled directly by the non-Abelian 
gauge field interactions in the fundamental forces of nature. Heavy-ion collisions at top RHIC and 
LHC energies produce strongly coupled plasma with a low value of fiR where lattice QCD predicts a 
smooth crossover between the QGP liquid and a hadron resonance gas. Experimentally mapping the 
QCD phase diagram is one of the big unsolved challenges in the field. A first Beam Energy Scan (BES-I), 
with Au+Au collisions at center-of-mass energies between 39 and 7.7 GeV to explore QCD matter at 
baryon chemical potentials 110 MeV < fj>B < 420 MeV, was completed in 2014. BES-I led to a number 
of intriguing observations of non-monotonic beam energy dependences of several flow and fluctuation 
observables which might be connected with the appearance of a first-order phase transition at large hb- 
State-of-the-art lattice QCD calculations combined with dynamical modeling, using a hybrid approach 
that couples viscous fluid dynamics for the QGP liquid with a microscopic approach to the critical phase 
transition dynamics and the subsequent evolution of the hadronic phase, will be required to test these 
interpretations. A second Beam Energy Scan (BES-II) planned for 2018-2019, with significantly improved 
beam luminosity and upgraded detector capabilities, and concurrent improvements of the theoretical 
modeling at lower beam energies are needed to solidify the suggestive results from BES-I with precision 
measurements in the targeted energy region identified in BES-I. Unambiguous discovery of a critical point 
in BES-II would warrant additional measurements at a later time to further quantify its properties. 

Forward rapidity studies at high energy and the Color Glass Condensate 

Both RHIC and the LHC are capable of probing new, unmeasured physics phenomena at low longitudinal 
momentum fraction x. Data from the 2013 p+Pb run at the LHC will make it possible to study previously 
unreachable phase-space in the search for parton saturation effects (the Color Glass Condensate). Forward- 
rapidity detector upgrades in STAR and PHENIX will open the door to studies of saturation physics 
at RHIC. However, a complete exploration of parton dynamics at low x will require an Electron Ion 
Collider (E 1C). While a future EIC will deliver crucially missing precise information on the nuclear parton 
distribution functions in a kinematic regime where in heavy-ion collisions saturation effects are difficult 
to separate from QGP physics, forward rapidity studies in p+A and A+A collisions at RHIC and LHC 
provide access to low-a; physics in a complementary kinematic range. 

2.3 Recommendations 

2.3.1 New questions and opportunities in relativistic heavy-ion collisions 

Over the past decade, through a panoply of measurements made in heavy-ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), in concert with theoretical advances 
coming from calculations done using many different frameworks, we have obtained a broad and deep 
knowledge of what hot QCD matter does, but we still know little about how it works. These collisions 
create exploding little droplets of the hottest matter seen anywhere in the universe since it was a few 
microseconds old. We have increasingly quantitative empirical descriptions of the phenomena manifest in 
these explosions, and of some key material properties of the matter created in these “Little Bangs” which 
turns out to be a strongly coupled liquid. However, we still do not know the precise nature of the initial 
state from which this liquid forms, and know very little about how the properties of this liquid vary across 
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its phase diagram or how, at a microscopic level, the collective properties of this liquid emerge from the 
interactions among the individual quarks and gluons that we know must be visible if the liquid is probed 
with sufficiently high resolution. 

Answering these and other questions requires an intensive modeling and computational effort to simulta¬ 
neously determine the set of key parameters needed for a multi-scale characterization of the QGP medium 
and the initial state from which it emerges. This phenomenological effort requires broad experimental 
input from a diverse set of measurements, including 1) the completion of the heavy quark program to 
measure the diffusion coefficient of heavy quarks, 2) energy scans to map the phase diagram of QCD and 
the dependence of transport coefficients on the temperature and baryon number chemical potential, 3) 
collisions of nuclei with varied sizes, including p+A and very high multiplicity p+p collisions, to study the 
emergence of collective phenomena, 4) the quantitative characterization of the electromagnetic radiation 
emitted by the Little Bangs and its spectral anisotropies, and 5) a detailed investigation at RHIC and 
LHC of medium effects on the production rates and internal structure of jets of hadrons, for multi-scale 
tomographic studies of the medium. These considerations lead us to our 


Recommendation #1: 

The discoveries of the past decade have posed or sharpened questions that are central 
to understanding the nature, structure, and origin of the hottest liquid form of matter 
that the universe has ever seen. As our highest priority we recommend a program to 
complete the search for the critical point in the QCD phase diagram and to exploit the 
newly realized potential of exploring the QGP’s structure at multiple length scales with 
jets at RHIC and LHC energies. This requires 

• implementation of new capabilities of the RHIC facility (a state-of-the-art jet 
detector such as sPHENIX and luminosity upgrades for running at low energies) 
needed to complete its scientific mission, 

• continued strong U.S. participation in the LHC heavy-ion program, and 

• strong investment in a broad range of theoretical efforts employing various analytical 
and computational methods. 


RHIC and the LHC, together, provide an unprecedented opportunity to study the properties of QCD 
matter. While collisions at the LHC create temperatures well above those needed for the creation of 
QGP and may thus be able to explore the expected transition from a strongly coupled liquid to a weakly 
coupled gaseous phase at higher temperatures, the RHIC program enables unique research at temperatures 
close to the phase transition. Moreover, the unparalleled flexibility of RHIC makes possible collisions 
between a variety of different ion species over a broad range in energy. The combined programs permit 
a comprehensive exploration of the QCD phase diagram, together with precise studies of how initial 
conditions affect the creation and dynamical expansion of hot QCD matter and of the microscopic 
structure of the strongly coupled QGP liquid. There is no single facility in the short- or long-term future 
that could come close to duplicating what RHIC and the LHC, operating in concert, will teach us about 
Nature. 

2.3.2 A precision femtoscope to study the glue that binds us all 

It is generally believed that, at high energies, interacting gluon fields from the colliding hadrons dominate 
the energy deposition in the collision zone and the subsequent thermalization processes that eventually 
lead to the creation of a quark-gluon plasma. The possibility to precisely measure the gluon wave 
functions of the incoming nuclei, in order to complement the empirical determination of the gluonic initial 
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state from flow fluctuations and correlations in heavy-ion collisions, is expected to further solidify the 
determination of the key parameters characterizing the thermodynamic and transport properties of the 
QGP. The construction of an Electron Ion Collider (EIC) will address this need. 

The EIC will image the gluons and sea quarks in the proton and nuclei with unprecedented precision and 
probe their many-body correlations in detail, providing access to novel emergent phenomena in QCD. It 
will definitively resolve the proton’s internal structure, including its spin, and explore a new QCD frontier 
of ultra-dense gluon fields in nuclei at high energy. These advances are made possible by the EIC’s unique 
capability to collide polarized electrons with polarized protons and light ions at unprecedented luminosity 
and electrons with heavy nuclei at high energy. An EIC will be absolutely essential to maintain U.S. 
leadership in fundamental nuclear physics research in the coming decades: 


Recommendation #2: 

A high luminosity, high-energy polarized Electron Ion Collider (EIC) is the U.S. QCD 
community’s highest priority for future construction after FRIB. 


2.3.3 New initiatives to further strengthen nuclear theory 

Due to the complexity of heavy-ion collision dynamics, the success of the research program outlined in 
this document hinges on continued strong support of a broad range of theoretical activities. In addition 
to a healthy base program, the field requires specific support in high-performance computing and of 
collaborative efforts focused on addressing complex issues that need contributions from several scientists 
or groups of scientists with complementary expertise. Examples for the first need are the determination 
of thermal equilibrium properties and response functions of the QGP from lattice QCD, the dynamical 
simulation of the thermalization processes in the initial stage of heavy-ion collisions, and comprehensive 
model/data comparisons involving large data sets and complex dynamical models that aim at extracting 
key model parameters with quantified uncertainties. An example for the latter is the DOE-funded JET 
Topical Collaboration that addresses the problem of turning measurements of the medium modification of 
jets and electromagnetic probes into precise tomographic tools that yield quantitative information about 
the properties and dynamical evolution of the dense QCD matter created in heavy-ion collisions. 

The high-performance computing needs of the nuclear theory community were addressed at a Town 
Meeting in Washington, D.C., on July 14-15, 2014, which passed the following recommendation and 
request: 

Recommendation: Realizing the scientific potential of current and future experiments demands large-scale 
computations in nuclear theory that exploit the US leadership in high-performance computing. Capitalizing 
on the pre-exascale systems of 2017 and beyond requires significant new investments in people, advanced 
software, and complementary capacity computing directed toward nuclear theory. 

Request: To this end, we ask the Long-Range Plan to endorse the creation of an NSAC subcommittee to 
develop a strategic plan for a diverse program of new investments in computational nuclear theory. We 
expect this program to include: 

• new investments in SciDAC and complementary efforts needed to maximize the impact of the 
experimental program; 

• development of a multi-disciplinary workforce in computational nuclear theory; 

• deployment of the necessary capacity computing to fully exploit the nations leadership-class 
computers. 

At the QCD Town Meeting at Temple University, the “Phases of QCD Matter" subcommunity endorsed 
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this resolution and amended it as follows: 


Recommendation #3: 

We endorse the new initiatives and investments proposed in the Recommendation and 
Request received from the Computational Nuclear Physics Town Meeting, at a level to 
be determined by the requested NSAC subcommittee. In addition, we recommend new 
funding to expand the successful “Topical Collaborations in Nuclear Theory” program 
initiated in the last Long Range Plan of 2007, to a level of at least one new Topical 
Collaboration per year. 


2.3.4 Educating and mentoring the next generation of scientists 

A continuous pipeline of highly and broadly educated young scientists is not only the lifeblood of the 
U.S. Nuclear Science program, but also a guarantor of our nation’s continued technological and economic 
strength and its ability to innovate. At a separate Town Meeting dedicated to Education and Innovation 
(Michigan State University, August 6-8, 2014) the following resolutions were passed: 

1. Education and mentoring of the next generation nuclear scientists as well as dissemination of 
research results to a broad audience should be recognized by all researchers as an integral part of 
the scientific enterprise. 

2. Nuclear science is an active and vibrant field with wide applicability to many societal issues. It is 
critical for the future of the field that the whole community embraces and increases its promotion 
of nuclear science to students at all stages in their career as well as to the general public. 

3. Researchers in nuclear physics and nuclear chemistry have been innovative leaders in the full 
spectrum of activities that serve to educate nuclear scientists as well as other scientists and the 
general public in becoming informed of the importance of nuclear science. The researchers are 
encouraged to build on these strengths to address some of the challenges in educating an inclusive 
community of scientists as well as those on the path to future leadership in nuclear science. 

4. The interface between basic research in nuclear physics and exciting innovations in applied nuclear 
science is a particularly vital component that has driven economic development, increased national 
competitiveness, and attracts students into the field. It is critical that federal funding agencies 
provide and coordinate funding opportunities for innovative ideas for potential future applications. 

In separate sessions at the QCD Town Meeting at Temple University, the "Phases of QCD Matter” and 
“QCD and Hadron Structure” subcommunities voted unanimously for the following 


Recommendation #4: 

The QCD community endorses and supports the conclusions from the Education and 
Innovation Town Meeting. 
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3 Quantifying the Properties of QCD Matter - 
Present Status and Future Opportunities 


3.1 Completing the Little Bang Standard Model 


One of the most important recent discoveries in heavy-ion collisions is that density fluctuations from the 
initial state of heavy-ion collisions survive through the expansion of the fireball showing up as correlations 
between produced particles |[3j-ll|. Prior work had mostly approximated the incoming nuclei as smooth 
spheres and the initial overlap region as an ellipse. The survival of density and geometry fluctuations 
was first hinted at in measurements of cumulants related to shape of the V 2 distribution [12,[l31. The 
picture started to become more clear after measurements were made in Cu+Cu collisions where the 
fluctuations were more prominent in the smaller system |14|. Ultimately, a new paradigm emerged as 
the structure of the initial state was found to play a central role in two-particle correlation functions 
and the previous measurements of V 2 were generalized to v n , a spectrum carrying information about 
both the initial densities in the collision as well as the dissipative properties of the subsequent plasma 
phase 1151. The survival of the initial state fluctuations is intimately related to the earlier finding that the 
QGP discovered at RHIC is the most perfect fluid known fl6f|l8|, with a viscosity to entropy density ratio 
77/s near the string theory bound 1191. The low viscosity plasma phase acts as a lens (albeit of strongly 
non-linear character), faithfully transferring the geometric structure of the initial density distributions, with 
its associated distribution of pressure gradients which act as a hydrodynamic force, into the final state. 
There it shows up most prominently as correlations between produced particles. Quantum fluctuations in 
the initial state cause these correlations to fluctuate from event to event. 


A standard dynamical framework for heavy-ion collisions 


Descriptions of these new phenomena have required the development of a new dynamical framework for 
heavy-ion collisions. It includes i) modeling of initial-state quantum fluctuations of nucleon positions and 
sub-nucleonic color charges and the color fields generated by them, ii) a description of the pre-equilibrium 
dynamics that evolves the initial energy-momentum tensor by solving either the (2+l)-dimensional 
Yang-Mills equations for the gluon fields (weakly-coupled approach) or Einstein’s equations of motion 
in five-dimensional anti-deSitter space (strongly-coupled approach), followed by iii) the rapid transition, 
event-by-event, to second-order viscous relativistic fluid dynamics, and iv) a late-stage hadron phase 
described by microscopic transport calculations. While there is widespread agreement on the general 
structure of such a standardized dynamical approach, it has not yet reached the level of uniqueness that 
would justify calling it the “Little Bang Standard Model” 1201. Model comparisons with experimental data 


that illustrate the state of the art in dynamical modeling can be found in 118 21 -291. With the existence 
of a reliable equation of state from lattice QCD calculations 30-33| a crucial degree of uncertainty in 
hydrodynamic modeling could be eliminated, enabling the development of a complete hydrodynamic 
space-time model. With this full space-time picture in hand, the comparisons of model calculations to 
harmonic decompositions of correlation functions (\A>f) at RHIC and the LHC (shown in Fig. [l| have 
reduced the uncertainty on 77/s by a factor of 10 1 34 1 . With this newfound precision, studies suggest that 
77/s is smaller for RHIC collisions (right panel of Fig. [I]) than it is at the LHC (left panel), consistent 
with a temperature dependent 77/s with a minimum near the critical temperature. In the next phase of 
study we seek to 1) accurately determine the temperature dependence of 77/s (aided by the Beam Energy 
Scan Program at RHIC) and 2) develop a clearer picture of the high density gluon fields that form the 
precursor of the plasma phase (aided by the p+A program and ultimately by an Electron Ion Collider). 

The path towards a Standard Model for the Little Bang 


So what is needed to turn this standard dynamical framework into the "Little Bang Standard Model”? 
One fundamental challenge along the way is the need to determine simultaneously the space-time picture 
of the collective expansion and the medium properties that drive this expansion 1201. A unique and reliable 
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Figure 1: Model calculations compared to measurements of the harmonic decomposition of azimuthal 
correlations produced in heavy-ion collisions. The left panel shows model calculations and data for v n vs. 
collision centrality in Pb+Pb collisions at yfsNN = 2.76 TeV. The right panel shows similar studies for 
the pt dependence of v n in 200 GeV Au+Au collisions. The comparison of the two energies provides 
insight on the temperature dependence of p/s. 


determination of these two unknowns is aided by measurements of multiple flow observables sensitive 
to medium properties in different stages of the evolution fl5 35 361. Due to the large event-by-event 


fluctuations in the initial state collision geometry, in each collision the created matter follows a different 
collective expansion with its own set of flow harmonics (magnitude v n and phases <J> n ). Experimental 
observables describing harmonic flow can be generally given by the joint probability distribution of the 
magnitudes v n and phases 4> n of flow harmonics: 
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dN f 


evts 


-^evts dVndVjyi . . . d*& n d*& r 
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Specific examples include the probability distribution of individual harmonics p(v n ), flow de-correlation 
in transverse and longitudinal directions, and correlations of amplitudes or phases between different 
harmonics {p(v n ,v m ) or p(4> n , 4> m )). The latter are best accessed through measurements of correlations 
with three or more particles. The joint probability distribution ([!]) can be fully characterized experimentally 
by measuring the complete set of moments recently identified in Ref. 1371. With the added detail provided 
by these measurements, hydrodynamic models can be fine-tuned and over-constrained, thus refining 
our understanding of the space-time picture and medium properties of the heavy-ion collisions. Initial 
measurements of some of these observables 


and comparison to hydrodynamic models 125 41 -431 
already provided unprecedented insights on the nature of the initial density fluctuations and dynamics 
of the collective evolution. Flowever, at this point none of the state-of-the-art hydrodynamic models 


properly accounts for the dynamical fluctuations generated during the evolution by thermal noise 44-47 
- future quantitative work will need to address these, too. 


Precision determination of key QGP parameters 

The agreement between the model and the data shown in Figure [T| suggests that the essential features 
of the dynamic evolution of heavy-ion collisions are well described by current models. These model 
calculations depend on a significant number of parameters that are presently poorly constrained by 
fundamental theory, and a reliable determination of the QGP properties requires a systematic exploration 
of the full parameter space. An example of such an exploration |48,491 is shown in Figure [2] where the 
shape of the QCD EOS is treated as a free parameter. The left panel shows a random sample of the 
thousands of possible Equations of State, constrained only by results on the velocity of sound obtained by 
perturbative QCD at asymptotically high temperature and by lattice QCD at the crossover transition 
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Constraining Eg. of State with RHIC/LHC Data (MADAI Collab.) 
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Figure 2: The QCD Equation of State (more precisely: the squared speed of sound, c 2 , as a function of 
temperature) from lattice QCD calculations and from models constrained by data from RHIC and the 
LHC. The red lines delimit the present uncertainty range of c 2 from lattice QCD. The blue lines show a 
number of parametrizations of c 2 (T)\ in the left panel, they are constrained only by their asymptotic 
values at T = 165MeV ~ T c and T = oo, in the right panel they are additionally required to provide an 
acceptable fit to experimental data as described in |48,491. The right panel shows that experimental data 
give preference to an Equation of State consistent with lattice QCD. This demonstrates that our model 
of the collision dynamics is good enough to allow us to study emergent properties in QCD matter. 


temperature. They are compared to the EOS determined from lattice QCD 1331.The right panel shows a 
sample of the Equations of State allowed by experimental data. The results of this study suggest that 
data at RHIC and the LHC require an EOS consistent with that expected from QCD. This demonstrates 
that our model of heavy-ion collisions describes the dynamics of the collisions well enough that we can 
extract information on the emergent properties of finite temperature QCD from the experimental traces 
left by the tiny droplet of QGP created in the collisions. These state-of-the-art models can therefore be 
used to both determine properties of finite temperature QCD currently inaccessible to lattice calculations 
and to provide an accurate space-time profile needed for modeling other processes like jet quenching. 
Figure [3] shows a schematic representation of our current uncertainty on the temperature dependence of 
r//s in QCD matter. While many of the existing measurements are accurate enough, as seen in Fig. [l] 
to determine r//s with much greater precision if all other model parameters were already known, the 
non-linear simultaneous dependence of the observables on multiple parameters does not yet allow one to 
translate the high quality of these experimental data into a more precise estimate of 77 /s. The studies 
shown in Figures [2] and [l]suggest, however, that a more complete set of measurements of the moments 
of the joint probability distribution ([Tj) at the LHC and RHIC (particularly in the Beam Energy Scan), 
coupled with extensive quantitative modeling, will provide the desired access to ( rj/s)(T ) in and around 
the transition temperature where hadrons melt into a Quark Gluon Plasma, and strongly reduce the width 
of the blue uncertainty band in Fig. [3| 

The temperature dependence of 77 /s can be further constrained by measuring directly emitted photons 
and dileptons. Their lack of strong interactions makes them penetrating probes that reflect the medium 
properties at the time of their emission which, on average, precedes that of the much more abundant 
strongly interacting hadrons. The slopes of their transverse momentum spectra and their anisotropic flow 
coefficients therefore show sensitivity to the hydrodynamic flow and dissipative medium properties at 
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higher average temperature than the corresponding hadronic observables 150,511. Recent elliptic and 
triangular flow measurements of direct photons in Au+Au collisions at RHIC 1521 and Pb+Pb collisions 
at the LHC 1531 have presented a puzzle to theorists because the measured flow anisotropies are much 


larger than theoretically predicted 150 511. This may point to stronger electromagnetic radiation from 
the critical crossover region between the QGP and hadronic phases 1541 than presently thought. Future 
higher-luminosity runs and improved detector capabilities will provide much more precise data for these 
very difficult measurements. In concert with state-of-the-art dynamical modeling, this is expected to not 
only yield a resolution of this “direct photon flow puzzle” but also lead to tighter constraints on the 
transport properties of QCD matter at higher temperatures. 



T[GeV] 


Figure 3: The temperature depen¬ 
dence of the specific viscosity rj/s. 
The blue band represents the range 
allowed by our current understand¬ 
ing based on comparing experimen¬ 
tal data to models with a mini¬ 
mum at the transition temperature. 
pQCD calculations and the string 
theory lower bound are also shown. 
The shaded vertical regions repre¬ 
sent the ranges of initial tempera¬ 
tures probed by RHIC and the LHC. 


The origin of collective behavior: Understanding the initial state and its thermalization 

Completion of the Little Bang Standard Model will also require a better understanding of the gluon 
densities in the initial overlap regions. As discussed in sections 3.4 and [4] this topic can be studied in p+A 
collisions and ultimately at an electron ion collider (EIC). The physics of p+A collisions has proven to be 
richer than originally anticipated. While collective flow phenomena are firmly established in sufficiently 
central collisions of heavy nuclei, these effects were not expected in p+p and p+A collisions. It was 
widely assumed that, as the mean-free-path of the matter approaches the characteristic size of the system, 
the effects of viscous damping would become too strong and invalidate a hydrodynamic description. 
Surprisingly, correlations that are long-range in rapidity and similar to those measured in A+A collisions 
have now also been observed at the LHC in rare high-multiplicity p+p collisions 1551 (corresponding to 
high gluon-density initial states). Subsequent measurements revealed similar phenomena at both the LHC 
and RHIC in high multiplicity p+Pb 


56-58 1 and d+Au 1591 collisions. In particular, the observed mass 
ordering of the spectral slopes 160,61 and of v n (pT ) |62)-|67| is reminiscent of the effects arising from a 
common radial flow boost in A+A collisions. Multi-particle measurements 168,69 1 show unambiguously 
that the correlations in high-multiplicity p+Pb collisions are collective. 


Understanding the relationship between initial and final-state effects in these systems appears to be 
non-trivial and has now developed into a very active research area. While hydrodynamic models with 
strong final-state interactions may provide a natural interpretation for many of the observed features, 
their applicability in such small systems, including the process of rapid thermalization, requires additional 
scrutiny 1701. Meanwhile, other novel mechanisms, exploiting quark and gluon momentum correlations in 
the initial state, have been proposed as alternative interpretations of the observed long-range correlations 
and even predicted the correlation results in p+p collisions. Recently collected data from ! He+Au 
collisions will shed additional light on this question, as will improved correlation measurements at forward 
rapidity, enabled by current and future detector upgrades. These programs play an important role in 
1) understanding the initial-state component of our Little Bang Standard Model and its subsequent 
modification by final-state interactions, 2) testing how small a system can be and still exhibit fluid-like 
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phenomena, 3) providing new opportunities to probe the spectrum of fluctuations in high gluon density 
matter, and 4) mapping the transition to a classical description of gluonic matter at high density. 

It is suspected that the rapid formation of almost perfectly liquid hot QCD matter in heavy-ion collisions 
may be related to the emergence of universal characteristics in high-density gluon matter at zero 
temperature that is predicted to dominate the low-.x component of the nuclear wave function when probed 
at high energy. To explore this connection, precision measurements of the nuclear wave function at an 
EIC will be required to complement nuclear collision experiments with small and large nuclei. 


3.2 Mapping the QCD phase diagram 


When the first protons and neutrons and pions formed in the microseconds-old universe, and when they 
form in heavy-ion collisions at the highest RHIC energies and at the LHC, they condense out of liquid 
quark-gluon plasma consisting of almost as much antimatter as matter. Lattice calculations 171—-731 show 
that QCD predicts that, in such an environment, this condensation occurs smoothly as a function of 
decreasing temperature, with many thermodynamic properties changing dramatically but continuously 
within a narrow temperature range around the transition temperature T c 6 [145 MeV, 163 MeV] 133 , 731 , 
referred to as the crossover region of the phase diagram of QCD, see Fig. [4j In contrast, quark-gluon 



Figure 4: A sketch illustrating the experi¬ 
mental and theoretical exploration of the 
QCD phase diagram. Although experi¬ 
ments at highest energies and smallest 
baryon chemical potential are known to 
cross from a QGP phase to a hadron gas 
phase through a smooth crossover, lower 
energy collisions can access higher baryon 
chemical potentials where a first order 
phase transition line is thought to exist. 


plasma doped with a sufficient excess of quarks over anti-quarks may instead experience a sharp first 
order phase transition as it cools, with bubbles of quark-gluon plasma and bubbles of hadrons coexisting 
at a well-defined critical temperature, much as bubbles of steam and liquid water coexist in a boiling 
pot. The point where the doping of matter over antimatter (parametrized by the net baryon number 
chemical potential hb ) becomes large enough to instigate a first order phase transition is referred to 
as the QCD critical point. It is not yet known whether QCD has a critical point 174-781, nor where 
in its phase diagram it might lie. Lattice calculations become more difficult or more indirect or both 
with increasing hb and, although new methods introduced within the past decade have provided some 
hints 175 77 791, at present only experimental measurements can answer these questions definitively. The 


theoretical calculations are advancing, however, with new methods and advances in computational power 
both anticipated. 

The phase diagram of QCD, with our current knowledge schematically shown in Fig. [4| is the only 
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phase diagram of any form of matter in Nature that we have a chance of both mapping experimentally 
and relating directly and quantitatively to our fundamental description of Nature, the Standard Model. 
With QCD the only strongly interacting theory in the Standard Model, mapping the transition region 
of its phase diagram is a scientific goal of the highest order. In the long term, successfully connecting 
a quantitative, empirical understanding of its phases and the transitions between phases to theoretical 
predictions obtained from the QCD Lagrangian could have ramifications in how we understand phases of 
strongly coupled matter in many other contexts. 

RHIC’s unique capability to probe the QCD phase diagram 

A major effort to use heavy-ion collisions at RHIC to survey the phase diagram of QCD is now underway. 
The excess of matter over antimatter in the exploding droplet produced in a heavy-ion collision can be 
increased by decreasing the collision energy, which reduces the production of matter-antimatter symmetric 
quark-antiquark pairs and gluons relative to the quarks brought in by the colliding nuclei, thus increasing 
[is- Decreasing the collision energy also decreases the maximum, i.e. initial, temperature reached by 
the matter produced in the collision. A series of heavy-ion collision measurements scanning the collision 
energy 1801 can therefore explore the properties of matter in the crossover region of the phase diagram, 
matter that is neither quark-gluon plasma nor hadronic or both at the same time, as a function of the 
doping {is- Such a program can scan the transition region of the QCD phase diagram out to fi b values 
that correspond to collision energies below which the initial temperature no longer reaches the transition. 
If the crossover region narrows to a critical point within this experimentally accessible domain, an energy 
scan can find it. RHIC completed the first phase of such an energy scan in 2014, taking data at a series 
of energies (^snn = 200, 62.4, 39, 27, 19.6, 14.5, 11.5 and 7.7 GeV) corresponding to values of hb that 
range from 20 to 400 MeV. Data from these experiments at RHIC 180,811 and from previous experiments 
confirm that lower-energy collisions produce QGP with higher / ib , as anticipated. A selection of data 
from BES-I that exhibit interesting non-monotonic behavior as a function of collision energy is shown in 
Hg.0 

RHIC is, and will remain, the optimal facility in the world for studying matter in the crossover region 
and searching for a possible critical point in the so far less well understood regions of the phase diagram 
with larger fj, b■ What makes RHIC unique is both its wide reach in hb and that it is a collider, meaning 
that the acceptance of detectors, and hence the systematics of making measurements, change little as a 
function of collision energy. Accelerator and detector performance has been outstanding during the first 
phase of this program, referred to as Beam Energy Scan I or BES-I. Measurements of all the important 
observables targeted in the planning of this campaign have now been made in collisions with energies 
varying by a factor of 25, allowing for a first look at a large region of the phase diagram of QCD. 

The need for BES-II and accompanying advances in theory 


We still await with interest the results from the most recent run in this program at ^/snn = 14.5 GeV, where 
data were taken only a few months ago; and for a number of important observables the measurements 
made at and below ^/snn = 19.6 GeV have quite limited statistics. Nevertheless, it is already possible 
to see trends and features in the data that provide compelling motivation not only for a strong and 
concerted theoretical response aiming at quantitative precision, but also for much higher statistics data 
at the lower end of the energy range (i.e. at the largest achievable values of ib ) that will be provided by 
the second phase of the Beam Energy Scan program (BES-II) in 2018 and 2019. The goals of BES-II, 
as described in more detail below, are to follow through and turn trends and features into definitive 
scientific conclusions. To this end, the accelerator physicists at RHIC are planning a machine upgrade 


to provide electron cooling to increase the beam luminosity at these energies by about a factor 10 1801 


Targeted new detector capabilities (an Inner Sector Upgrade for the STAR TPC and a new event plane 
detector) will also increase the sensitivity to signals already observed in BES-I and open the door to new 
measurements described below in the BES-II campaign 


Experimental discovery of a first order phase transition or a critical point on the QCD phase diagram 
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Figure 5: Three selected observables that all show 
interesting non-monotonic behavior as functions 
of collision energy around ^/snn ~ 15—20 GeV. 
Top panel: Rl ut —R 2 si de' measured via two-pion 
interferometry by STAR 1821, PFIENIX 1831, and 


ALICE 1841, reflects the lifetime of the collision 


fireball. Ri u t~Rside was predicted 85 to reach a 
maximum for collisions in which the hydrodynamic 
fluid forms at temperatures where the equation of 
state is softest. 

Middle panel: The rapidity-slope of the net pro¬ 
ton directed flow v\, dv\/dy. This quantity is 
sensitive to early pressure fields in the medium. 
Bottom panel: The kurtosis of the event-by¬ 
event distribution of the net proton (i.e. proton 
minus antiproton) number per unit of rapidity, 
normalized such that Poisson fluctuations give a 
value of 1. In central collisions, published results 
in a limited kinematic range 1861 show a drop be¬ 
low the Poisson baseline around ^/snn =27 and 
19.6 GeV. New preliminary data over a larger px 
range [87], although at present with substantial 
uncertainties, hint that the normalized kurtosis 
may, in fact, rise above 1 at lower yfsNN, as ex¬ 
pected from critical fluctuations 1881. The grey 
band shows the much reduced uncertainties an¬ 
ticipated from BES-II in 2018-2019, for the 0-5% 
most central collisions. 


would be a landmark achievement. The first goals of the BES program, however, relate to obtaining a 
quantitative understanding of the properties of matter in the crossover region of the phase diagram as it 
changes with increasing ps- Available tools developed over the last few years now make a quantitative 
comparison between theory and experiment tractable in the ^^-range below any QCD critical point. 
Success in this, in and of itself, would constitute a major and lasting impact of the RHIC program. 
Questions that can be addressed in this regime include quantitative study of the onset of various signatures 
associated with the presence of quark-gluon plasma and of the onset of chiral symmetry restoration as one 
traverses the crossover region. Data now in hand from BES-I provide key inputs and impetus toward this 
goal. Flere we give four examples, intended to be illustrative, of areas where a coherent experimental 
and theoretical effort is expected to have substantial impact on our understanding of QCD. In each case 
we also note the substantial impact expected from the additional measurements anticipated during the 
BES-II: 

1. The directed flow observable dvi/dy for net protons has been found to feature a dip as a function of 
collision energy (see middle panel in Fig. [5]), with a minimum at energies somewhere between ^/snn = 11.5 
and 19.6 GeV 189]. This has long been predicted in qualitative terms as a consequence of the softening of 
the equation of state in the transition region of the phase diagram 190,911. Several theoretical groups 
around the world have now begun hydrodynamic calculations with nonzero baryon density, deploying all the 
sophistication that has been developed very recently in the analysis of higher energy collisions, including 
initial fluctuations and a hadronic afterburner, in applications to these lower energy collisions. These 
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hydrodynamic+hadronic cascade calculations will be used to compare the dvi/dy data with equations of 
state in the crossover region of the phase diagram obtained from lattice calculations via Taylor expansion 
in hb/T 1921. This is a program where a quantitative comparison, successful or not, will be of great 
interest, since failure to describe the data could signal the presence of a first order phase transition. 
The precision of a comparison like this will be substantially improved in 2018-19 when BES-II data will 
allow dv\/dy to be measured for the first time with tightly specified centrality; the statistics available in 
the BES-I data sets limit present measurements to averages over collisions with widely varying impact 
parameters |89|. 

2. A second goal of the hydrodynamic calculations referred to above will be to use identified particle 
BES-I V 2 data to map, in quantitative terms, where and how hydrodynamics starts to break down at 
lower collision energies, and where, to an increasing extent, V 2 develops during the hadron gas phase 
when viscosities are not small, i.e. where the contribution of the partonic phase to observed measures of 
collectivity decreases in importance. A key future experimental input to this program is the measurement 
of the elliptic flow V 2 of the (/>-meson, which will be obtained with substantially greater precision in the 
BES-II program. The first measurements of V 2 of Q baryons at these collision energies, also anticipated 
in BES-II, will represent a further, substantial advance. Seeing cp mesons flowing like lighter mesons and 
Q baryons flowing like lighter baryons in collisions at a given energy would indicate that the dominant 
contribution to the collective flow in those collisions was generated during the partonic phase 1931. 

This component of the BES program, together with the following one, will yield guidance as to what the 
lowest collision energies are at which temperatures in the transition region of the phase diagram can be 
explored. That is, they will tell us out to what ys it will be possible for heavy-ion collisions, anywhere, 
to study matter in the crossover region and search for a possible critical point. 


3. Heavy-ion collisions at top RHIC energies and at the LHC have now seen several experimental 
phenomena (94-96 that may be related to the chiral magnetic effect (CME 197,981, see Sec ]3.5.2] ). 
In each case, alternative explanations are also being considered 199,1001. One of the intriguing BES-I 
results is that the three-particle correlations that are related to charge separation across the reaction 
plane, possibly induced by the CME, are robustly observable over most of the BES range but then seem 
to turn off at yfsNN = 7.7 GeV [1011. where the elliptic flow V 2 is still robust. This is an indication 
that r^-induced backgrounds alone do not explain the observed correlations. The observation that these 
three-particle correlations disappear at the lowest energy could prove crucial to understanding their origin 
and how they are related to the formation of QGP. On the theoretical side, lattice QCD calculations 
probing the response of the equation of state and transition temperature to the presence of external 
magnetic fields 1102—1041, as well as hydrodynamic calculations incorporating magnetic fields and chiral 
effects, are needed and are being pursued by several groups. On the experimental side, higher statistics 
BES-II data will make it possible to determine with much greater precision the snn at which this 
effect turns off and will also make it possible to measure the (related but theoretically more robust) 
chiral magnetic wave phenomenon 1105,1061, which has also been seen at top RHIC energy and at the 
LHC 1107 108), and which should turn off at the same ^/snn if these interpretations are correct. 

4. Theoretical developments over the past decade have identified specific event-by-event fluctuation 
observables most likely to be enhanced in collisions that cool in the vicinity of the critical point 1109,1101. 
Higher moments of the event-by-event distribution of the number of protons, or the net proton number, 
are particularly sensitive 1110-1121. STAR has now measured the first four moments (mean, variance, 
skewness and kurtosis) of the event-by-event distribution of net proton number and net charge at the 
BES-I energies 186,1131. At the lowest collision energies, although the statistics are at present rather 
limiting, there are interesting trends, including e.g. the drop in the kurtosis of the net-proton distribution 
at yj snn = 27 and 19.6 GeV (see bottom panel in Fig. [5]). This drop in and of itself can be at least 
partially reproduced via prosaic effects captured in model calculations that do not include any critical 
point. Theoretical calculations of the contributions from critical fluctuations predict 1881 that if the 
freezeout \ib scans past a critical point as the beam energy is lowered, this kurtosis should first drop 
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below its Poisson baseline and then rise above it. Both the drop and the rise should be largest in central 
collisions in which the quark-gluon plasma droplet is largest and therefore cools most slowly, allowing 
more time for critical fluctuations to develop 11141. A recent and still preliminary analysis 1871 of the 
data at y / SNN = 11.5 and 7.7 GeV over a larger range in px than measured before 1861, also shown in the 
bottom panel of Fig. [5] shows intriguing hints of a rise in the net proton kurtosis in central collisions, but 
the uncertainties are at present too large to draw conclusions. If this kurtosis does rise at t/snn values 
below 19.6 GeV, this would be difficult to understand in conventional terms and thus would be suggestive 
of a contribution from the fluctuations near a critical point. Determining whether this is so requires the 
higher statistics that BES-II will provide, as illustrated by the grey band in the bottom panel of Fig. [5] 


The present data on moments of both the net proton number and the net charge at the higher BES-I 
energies are already very useful, as they can be compared to lattice calculations of the Taylor expansions 
(in pb/T) of the baryon number and charge susceptibilities 11151. First versions of this comparison 
have been reported recently and are being used to provide an independent determination of how the 
freeze-out values of pb and T change with collision energy 1116-1191. Flowever, looking ahead, theoretical 
calculations will need to faithfully account for the dynamical evolution of the medium formed in the 
collision for a full quantitative exploitation of the experimental data. For the higher statistics BES-II 
data on the net proton kurtosis, skewness, and other fluctuation observables at low collision energies 
to determine the location of the critical point on the phase diagram of QCD, if one is discovered, or 
to reliably exclude its existence within the experimentally accessible region of the phase diagram, a 
substantial theoretical effort will be needed that couples the sophisticated hydrodynamic calculations 
referred to above with a fluctuating chiral and dynamically evolving order parameter. 


As the following fifth example illustrates, BES-II will also open the door to measurements that were not 
yet accessible in the first phase of the BES program: 

5. Dileptons are unique penetrating probes to study the chiral properties of hot and dense matter QL20,121 1 . 
The dielectron invariant mass distributions measured in the BES-I (in data taken at ^/snn = 200, 62.4, 
39, 27 and 19.6 GeV) have shown that there is a significant enhancement of low mass dileptons below 1 
GeV relative to a hadronic cocktail 1122—1251. The data to date are consistent with model predictions 
of a p resonance that broadens through interactions with the hadronic medium 11201 and melts in a 
partonic medium. This picture captures the essential features of the low-mass excess at small chemical 
potentials 11251. Data at lower energies with higher statistics are crucial in order to test the predicted 
dependence of dilepton yields on baryon density and draw firm conclusions. The dilepton measurements 
at and below ^/snn = 19-6 GeV that BES-II will provide will advance our understanding of the chiral 
properties of QCD matter with significant net baryon density. In addition, the low-mass excess can serve 
as a chronometer for the fireball lifetimes as the QCD phase diagram is scanned 11261. At higher masses 
(above 1 GeV) chiral restoration can be probed via the mixing between vector and axial-vector currents in 
the hot and dense environment, and there the slope of the dilepton mass spectrum provides access to the 
early temperature in the system, unaffected by collective flow 11261. 

Each of these five examples makes it clear that in order to maximize the physics outcome from BES-I 
and BES-II, a coherent effort between experimentalists and theorists working on QCD at nonzero T and 
Pb is essential and must be organized and supported. Indeed, there has been considerable progress in 
lattice QCD recently on the calculation of various QCD susceptibilities 1127,1281 and the QCD equation 
of state in the regime where pb is nonzero but sufficiently small compared to 3T C [I29j[l30). These 
lattice calculations provide the necessary inputs for extending to nonzero ps the kind of sophisticated 
hydrodynamic calculations (including initial fluctuations and a late stage hadron cascade) that have been 
developed over the past few years. For some purposes, these calculations additionally require dynamical 
coupling to a fluctuating chiral order parameter. 


In concert, such developments will provide the critical tools for obtaining from BES-I and BES-II data 
answers to fundamental questions about the phases, the crossover, and perhaps the critical point and 
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first order transition, in the QCD phase diagram. Quantitative understanding of the properties of matter 
in the crossover region where QGP turns into hadrons will come first. If there is a critical point with 
/ib < 400 MeV, BES-II data together with the theoretical tools now being developed should yield 
quantitative evidence for its presence. The span in T and q ,b that the flexibility of RHIC makes accessible, 
along with the mentioned technical advantages of measuring fluctuation observables at a collider and 
recent and planned detector and facility upgrades (e.g. low energy electron cooling), put RHIC in a globally 
unique position to discover a critical point in the QCD phase diagram if Nature has put this landmark in 
the experimentally accessible region. Late in the decade, the FAIR facility at GSI will extend this search 
to even higher n b if its collision energies continue to produce matter at the requisite temperatures. 

3.3 Probing hot QCD matter at multiple length scales 

Hard probes, such as jets and bottomonia, provide the unique opportunity to test the microscopic structure 
of hot QCD matter at characteristic momentum and length scales. Below we describe the qualitative 
and quantitative insights into QGP properties these probes have already provided and outline a future 
program that aims to understand how the QGP properties arise from its microscopic nature. 


3.3.1 Jets as microscopic probes of QGP 


The properties of jets 1311 and their emergence from pQCD have been extensively studied in high-energy 
physics |132||133| . One of the earliest discoveries at RHIC was the phenomenon of jet quenching , observed 
as a suppression of high-py hadrons and di-hadron correlations in Au+Au collisions 1134 1351. This 
can be understood as a medium-modification of jet showers 11361 through an enhanced rate of gluon 
bremsstrahlung, resulting in a depletion of high momentum parton fragmentation products. 

The virtuality of a hard parton within a jet establishes the intrinsic scale at which it resolves fluctuations 
in the medium. As partons cascade down to lower virtualities, they probe the medium over a multitude of 
length scales. As long as this resolution scale is much larger than A qcd, the parton will be weakly coupled 
with the medium and pQCD can be applied to describe its propagation 1137—1461. Some partons within a 
jet may lose sufficient amounts of energy to encounter non-perturbatively strong coupling 1147—1501. 

Jet transport parameters from single-hadron suppression 

The medium-induced changes to the shower radiation pattern can be described as longitudinal drag/diffusion, 
transverse diffusion, and enhanced splitting of the propagating partons. The transport coefficients q [151| 
and e 11521 quantify the transverse diffusion and longitudinal drag, respectively. 

Since the last NP long-range plan, enormous progress has been achieved in the quantitative determination 
of q and e. This has been made possible by new high precision data from RHIC and LHC, as well as 
coordinated theory efforts, notably by the JET collaboration. One now obtains a very good description of 
the combined RHIC and LHC data on single hadron suppression, as shown in the left panel of Fig. [6j 
The extracted values and temperature dependence of the dimensionless ratio g/T 3 are plotted in the 
right panel of Fig. [6] Using identical initial states and hydro simulations, the JET collaboration has 
carried out a systematic analysis of a wide range of pQCD based energy loss schemes 11531. Compared 
to earlier studies |154|, where the extracted q varied by an order of magnitude, the values in the new 
analysis 11531 differ at most by a factor of 2. This makes it possible for the first time to discern the 
medium temperature dependence of g/T 3 and demonstrates that quantitative properties of the QGP can 
be extracted from jet modification data. 

Fully reconstructed jets 

Fully reconstructed jets have been a crucial tool in high energy physics for precision tests of pQCD 11311. 
Experimental progress now allows one to isolate and reconstruct the entire jet shower from the high- 
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Figure 6: A comparison of several different pQCD based energy loss energy loss schemes to the measured 
leading hadron suppression in central events at RHIC and LHC 11551, and the extracted scaled (i.e. 
dimensionless) transport coefficient q/T 3 along with its dependence on temperature 11531. 

multiplicity underlying event in heavy-ion collisions. Some of the first LHC heavy-ion results included 
the observation 1157,1581 of highly asymmetric dijet events 1159-1621, which provided a striking visual 
demonstration of the energy loss. Since these initial measurements, experimental control over the 
measured jet energies and the understanding of the role of underlying event fluctuations have improved 
substantially, resulting in precise measurements of jet suppression and the properties of quenched jets. 
Photon-jet correlations 1163-1651 have made it possible to select the parton kinematics and flavor before 
energy loss, providing direct evidence for the degradation of the parton energy as they traverse the 
medium 1166-1681. 

These measurements at the LHC, and pioneering studies at RHIC 11691, have contributed qualitatively 
new information QL70-1721 on jet modification by the medium. Studies of the cone-size dependence of jet 
spectra 1173], jet-hadron correlations (174| and the overall energy flow in dijet events demonstrate that 
the energy transported out of the hard core of the jet, as seen in hadron suppression, is not contained 
within typical jet cone radii jl75,1761. Rather, the energy is recovered in low momentum modes (few GeV 
or less) at large angles to the jet direction. The jet structure within the cone reveals only a moderate 
modification, undergoing a softening and broadening of the fragmentation pattern. 

Heavy-quark jets and the mass-dependence of parton energy loss 

The LHC also made new tests of the mass dependence of energy-loss possible using heavy-quark mea¬ 
surements. Data on D-meson production and non-prompt J/tj} from B-meson decays (for px up to 20 
GeV/c) exhibit the predicted mass hierarchy [177,1781, with the heaviest quarks losing significantly less 
energy than the lighter flavors. At jet energies of 80 GeV and above, i.e. much higher than the quark 
mass, tagged 6-quark jets show a similar suppression pattern as inclusive jets, as theoretically expected. 
At RHIC, pioneering studies of electrons from heavy flavor decays 1179 , of directly reconstructed D° 
mesons |180|, and of bottom/charm-separated non-photopic electrons |18l][l82j have shown considerable 
heavy-flavor suppression, with nuclear modification factors approaching those of light-flavor hadrons in 
the range 5G eV/c<pT < 10 GeV/c. 


Future plans - facilities, detectors and measurements 

Future jet-based studies built on the achievements at RHIC and LHC will address fundamental questions 
about the nature of QGP. These include precise measurements of QGP transport coefficients as a function 
of temperature, a detailed characterization of the QGP response to the parton energy loss and studies of 
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Figure 7: Shown on the left are jet virtuality evolution paths simultaneous with the QGP temperature 
evolution for central Au+Au and Pb+Pb collisions at RHIC and the LHC. The vertical axis shows the 
inverse length scale probes in the medium and the blue region where the medium is substantially altering 
the vacuum parton splitting 11561. This displays the complementarity of RHIC and LHC jet probes in 
surveying distinct regions of possible QGP excitations and temperature. The right panel shows the 
momentum regions covered by penetrating probe measurements at RHIC, enabled by high luminosity and 
the sPHENIX upgrade, and at the LHC with full Run II and Run III statistics. 


the modification the jets' angular and momentum structure as a function of angular and momentum scale. 
In combination, the goal of these studies is to understand the microscopic (or quasi-particle) nature of 
QGP at varying scales, and to understand how the macroscopic QGP liquid emerges from the underlying 
QCD degrees of freedom. A schematic sketch of the present and expected future resolving power for the 
structure of QCD matter at different temperatures in RHIC and LHC experiments is shown in Fig. El 11561. 


This program will be enabled by the evolution of the RHIC and LHC accelerator facilities, upgrades to 
the existing experiments and the construction of a state-of-the-art jet detector at RHIC, sPHENIX 11561. 
In parallel, the recently emerged experiment/theory collaborations will be strengthened and expanded, to 
fully utilize the increased precision and range of experimental observables. 


At BNL, three heavy-ion running campaigns are currently envisioned to complete the RHIC mission, as 
shown in Fig. [8] (left). A key goal of the 2014-16 period is to measure heavy flavor probes of the QGP 
using the newly installed silicon vertexing detectors in PHENIX and STAR. For 2018-19, Phase II of the 
RHIC Beam Energy Scan program (BES-II) is foreseen. In 2021-22, precision jet quenching and quarkonia 
measurements will be made possible by the installation of sPHENIX. 

The proposed sPHENIX detector, shown in the right panel of Fig. [8] would dramatically extend the 
range of jets measurable at RHIC, as well as provide precision spectroscopy of quarkonia, in particular 
the upsilon states. The program relies on a very high data acquisition bandwidth, which combined with 
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Figure 8: (Left) The timeline for future RHIC and LHC heavy-ion running. (Right) A cut-away view 
of the proposed sPHENIX upgrade to the PHENIX experiment, showing the inner tracking system, the 
electromagnetic calorimeter, the BaBar solenoid and the hadronic calorimeter. 


RHIC II luminosities enables the measurements of jet energies up to 70 GeV. The sPHENIX design takes 
advantage of various technological advances to minimize costs. sPHENIX also provides one route for a 
for a smooth evolution to a full-capability EIC detector. 

The LHC is now preparing for Run II, foreseen to include p+p, p+Pb and heavy-ion data taking from 
2015 to 2018. Run II will be followed by a shutdown from 2018 to 2020 (LS2) and Run III from 2020 to 
2023. For both the p+p and Pb+Pb data taking, the LHC upgrades during the current shutdown should 
make collisions possible at close to the design energy, i.e., about 5 TeV per nucleon pair for Pb+Pb. 
Collision rates are expected to exceed the design values by up to one order of magnitude, with a combined 
Run II and III goal to deliver about 10 nb _1 per experiment. In combination, the increased collision 
energy and luminosity will increase statistics for high px probes by about a factor of 200. 

To exploit this improved accelerator performance, ALICE, ATLAS and CMS are undergoing significant 
upgrades during the current shutdown and LS2. For ATLAS and CMS, these upgrades are mostly driven 
by the needs of the p+p program. Of particular importance for future heavy-ion running are upgrades to 
the inner tracker detectors, extending the reach of charged particle measurements, and upgrades to the 
trigger systems, allowing the experiments to fully sample the expected collision data for hard probes. 

ALICE is preparing for Run II with an expansion of the calorimetric coverage (EMCAL) which will make 
possible dijet studies and improved jet triggering. During LS2 the experiment’s data taking capabilties 
will be significantly enhanced with major upgrades to detector readout and data aquisition systems. In 
addition a new silicon inner tracker will be installed. While precision measurements of low px open heavy 
flavor are the main physics driver for these upgrades, they also benefit the intermediate and high px jet 
quenching program. 

These RHIC and LHC facility and experiment upgrades benefit jet physics studies in three major ways. 
The statistical precision and kinematic reach for commonly used jet physics observables is vastly increased. 
For single charged hadrons and reconstructed jets, the px reach will be extended by a factor of 2-3, up 
to 40 GeV for hadrons and 70 GeV for jets at RHIC and 300 GeV and 1 TeV, respectively, at the LHC. 
This will further improve the precision of the extraction of e.g. q and its temperature dependence from 
model comparisons. 

A second benefit is improved access to rare, highly specific observables, such as isolated photon + jet 
correlations at RHIC and LHC, as well as Z°+jet correlations at the LHC. For example, LHC experiments 
will record more than 5 x 10 5 photons with p '(, > 60 GeV, compared to about 3000 in the Run I data sets. 

Finally, the very high statistics jet samples enable the measurement of a new generation of jet shape 
observables. There is intense activity in developing generalized jet structure variables for the pp program 
at the LHC to maximize the efficiency of discovery measurements by improving quark/gluon discrimination 
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and the tagging of boosted objects. Heavy-ion studies of the modification of the jet momentum and 
angular structure through medium interactions will benefit greatly from these developments. 

Jet physics outlook 

In combination, these developments will enable a coherent RHIC and LHC physics program employing 
well-calibrated common observables. As shown in Fig. [7] future RHIC and LHC measurements will both 
overlap with and complement each other. This program will have three main components: 

First, a combined global analysis of RHIC and LHC data will lead to a precise determination of the 
values and temperature dependence of QGP transport coefficients. Both the RHIC and LHC final states 
represent an integral of jet-medium interactions over the evolution of both jet and medium from initial to 
final state. To disentangle the temperature dependence from this evolution it will be essential to deploy 
directly comparable observables (theoretically and experimentally) in different QGP temperature regimes, 
with particular focus on the fraction of their evolution spent in vicinity to the phase transition region that 
will be passed by all collisions. 

The second component relies on the increased systematic and statistical precision afforded by new probes 
(e.g. photon-jet) to identify the medium response to the modified jet radiation and further elucidate the 
liquid nature of the medium in its response to local perturbations. 

Third, precision measurements of the angular and momentum structure of jets can be used to characterize 
the microscopic structure of the QGP, using hard partons as probes in “Rutherford scattering” off effective 
QGP constituents or quasi-particles. Studies include potential modifications to the back-to-back jet 
scattering distributions, as well as modifications of the intra-jet angular structure. For the latter, the 
correlated angular and momentum evolution of the jet from the initial scattering to the final hadronic 
structure probes a wide range of scales, opening a window to interactions of jet and QGP constituents 
between vacuum-like and in-medium cascade regimes. To elucidate the physics of this intermediate 
window, a systematic variation of both the jet conditions and medium conditions and dynamics is 
necessary. Combining RHIC and LHC measurements will make it possible to control the initial density 
and temperature (in particular in respect to their vicinity to the critical temperature) and expansion 
dynamics of the system. The different energy regimes and tagging of particular initial states (photon+jet, 
6-tagged jets, multi-jet events) will enable the selection of different or common jet populations in relation 
to different medium conditions. Success in this long-term endeavor will require a global analysis of a 
diverse set of RHIC and LHC data in an improved, well controlled theoretical framework matched to the 
experimental observables. 


3.3.2 Probing the QGP with bottomonia 


A key probe of the strength and nature of the color interactions in the QGP are precision measurements 
of bottomonia states. The three lowest bound states of b-b (T Is, 2s, 3s) have a range of binding energies 
and length scales that allow one to understand the temperature dependence of plasma screening [183,184 
that leads to a characteristic sequential suppression pattern 11851 for in-medium quarkonium bound 
states. Compared to charmonia, bottomonia have the advantage that recombination in the QGP and 
at the phase boundary is expected to be smaller because of the smaller numbers of 6-quarks in the 
system, even at LHC energies, thus providing a more direct measure of color screening through the 
suppression mechanism 186 187 . On the other hand, in p+p collisions only 0.1% of b-b pairs form 


bottomonia (compared to ~ 1% for charm) 11881. Transport calculations suggest that bottomonium 
regeneration at both RHIC and LHC is relatively small 11891, which facilitates the interpretation of 
observed medium modifications of bottonium production yields in terms of color screening and dissociation 
effects |186,187||189|. 

First detailed bottomonium measurements at the LHC from CMS show the expected ordering by their 


25 























size and binding energies in the suppression of these states as a function of collision centrality 11901. At 
LHC energies regeneration may affect bottomonium production, making bottomonium measurements at 
RHIC an important benchmark measurement. First measurements at RHIC without resolving the three 
states have also been completed 1191,1921. Critical to the effort to understand the screening interactions 
as a function of temperature are simultaneous precision measurements at RHIC and the LHC (probing 
cooler and hotter plasmas) in p+p, p+A, and A+A, with sufficient statistics to map out the centrality 
and pt dependencies. Increased luminosities at the LHC in Run II and Run III will provide these data. At 
RHIC, additional statistics will be provided by the STAR MTD upgrade; ultimately, for matched precision 
at RHIC and the LHC, the sPHENIX upgrade is required. 


3.4 A polarized p+A program for initial-state and low-x phenomena 

Four dominant questions have been identified that can be explored in a polarized p+p and p+A program: 
1) What is the nature of the spin of the proton? 2) How do quarks and gluons hadronize into final-state 
particles? 3) How can we describe the multidimensional landscape of nucleons and nuclei? 4) What is the 
nature of the initial state in nuclear collisions? The first question is discussed in the “cold QCD” Town 
Meeting summary and the RHIC-Spin White Paper 11931. The remaining questions share particularly 
strong synergies with the goals of the hot QCD community and are discussed below. 

At the high collision energies of RHIC and the LHC, the available phase space for radiating small-x 
gluons and quark-antiquark pairs is very large. Since each emitted parton is itself a source of additional 
gluons, an exponentially growing cascade of radiation is created which would eventually violate unitarity. 
However, when the density of partons in the transverse plane becomes large enough, softer partons 
begin to recombine into harder ones and the gluon density saturates. This limits the growth of the 
cascade and preserves the unitarity of the S-matrix. The transverse momentum scale below which these 
nonlinear interactions dominate is known as the saturation scale Q s . The saturation scale grows with 
energy, but also with the size of the nucleus as Q ^ ~ A 1 / 3 . For high enough energies Q s is large and the 
corresponding QCD coupling is weak: a s (Q s ) <C 1. This makes it possible to calculate even bulk particle 
production using weak coupling methods, although the calculation is still nonperturbative due to the large 
field strengths. Because the gluonic states have large occupation numbers, the fields behave classically. 
The classical field theory describing the saturation domain is known as the “Color Glass Condensate” 
(CGC) 11941. 

The ideal probe of the CGC description are dilute-dense collisions, where a simple small projectile collides 
with a large nucleus. At RHIC and the LHC this makes p+A collisions important for understanding 
saturation. Significant progress has been made in describing the systematics of particle production as a 
function of transverse momentum and rapidity in p+p and p+A collisions with CGC calculations, which 
match the collinearly factorized perturbative QCD description at high transverse momenta. The case 
of saturation effects in multiparticle correlations as a function of azimuthal angle and rapidity remains 


more open (for more discussion of the ridge in p+p and p+A collisions see Sec. 3.1). While there are 
contributions to these correlations that originate already in the nuclear wave functions 11951, experimental 
evidence points to strong collective behavior also in the final states of p+A and even p+p collisions. The 
versatility of RHIC to systematically change the size of the projectile nucleus and complement p+A with 
d+A, 3 He+A etc. collisions over a wide range of collision energies is unparalleled and a key to exploring 
where these collective effects turn on. Precise and controlled access to the high energy nucleus can be 
used to disentangle the effects of strong collectivity in the initial wave functions and the final state. 

As noted elsewhere, an EIC |2| which can measure the transverse and longitudinal structure of the small-x 
gluons in nuclei is crucial for a deep understanding of hadron wave functions. However, measurements 
of forward photon, J /ip, Drell-Yan, inclusive and di-jet, and hadron/jet correlation probes in p+A 
collisions 1196, 197] provide a unique opportunity to make timely progress in our understanding and 
to complement the eventual EIC measurements. These measurements take advantage of the unique 
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capabilities afforded by the RHIC accelerator and detector complex, complemented with detector upgrades 
including the forward calorimetric system and forward tracking system required to carry out the proposed 
physics program at forward rapidities. In particular, the observables mentioned are related to different 
transverse momentum dependent gluon distributions (TMDs) that have previously been studied at higher 
x. Recent theoretical advances 11981 have clarified the relation of these distributions to each other via 
the CGC picture, in both p+A and e+A collisions, and they are opening up new exciting connections 
to polarized observables. One particular example of these connections is the possibility of extracting 
the saturation scale Q s in nuclei by comparing transverse single spin asymmetries measured in polarized 
p+p and polarized p+A collisions, for different nuclei and at different beam energies 1199,2001. Another 
exciting opportunity in polarized p+A collisions is the possibility of extracting first information about 
generalized parton distributions (GPDs) for gluons. This can be achieved for instance in exclusive J/T' 
production via photon-gluon fusion of quasi-real ultraperipheral photons from the nucleus with gluons 
from the polarized proton 11931. 


The study of these novel TMDs and GPDs will deepen our understanding of the momentum and spatial 
structure of polarized and unpolarized quarks and gluons in hadrons. These studies require as key 
ingredients establishing factorization and universality that help separate the structure of the hadron wave 
function from the dynamics of the probe. From this perspective, even if measurements at a future EIC 
with an electron probe provides unmatched precision, a polarized proton-nucleus program provides a 
complementarity that may prove essential. 



Figure 9: Color fields of the two nu¬ 
clei before the collision, from 12011. 


The theoretical description of the initial stage of quark gluon plasma formation has become increasingly 
detailed. State of the art calculations 125, 2021 now combine a fluctuating nuclear geometry with a 
microscopic QCD description of the dynamics of matter formation (see Fig. [9]), going beyond the Glauber 
model descriptions of the geometry. As discussed in Sec. 3T these initial-state calculations, combined 
with detailed measurements of correlations and fluctuations in the observed flow patterns have helped 
to significantly improve the precision of the first quantitative experimental determinations of e.g. the 
viscosity/entropy ratio rj/s. In this context, p+A collisions provide crucial consistency checks for our 
understanding of heavy-ion collisions: Clearly, theoretical descriptions developed for partonic interactions 
in a hot and dense QCD medium must also be consistent with the effects of the cold QCD medium on 
hadronization of a colored probe that can be studied in p+A collisions. Such collisions, across a range of 
beam energies and target species, will therefore provide important control experiments for our theoretical 
understanding of jet quenching in heavy-ion collisions, on a time scale that precedes the construction of 
an EIC. As noted previously, p+A collisions have provided surprises in their own right-we now understand 
the resolution of these to be sensitive to the detailed spatial structure of partons in both protons and 
heavier nuclei. The importance of a polarized p+A program is therefore two-fold: (i) It will provide 
unique and essential information on the parton structure of proton and nuclear wave functions, (ii) The 
implementation of this information in models of heavy-ion collision will provide more sensitive tests of 
and precise extraction of the parameters of the Little Bang Standard Model. 
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3.5 Opportunities and challenges in theory 

3.5.1 Progress since the 2007 Nuclear Physics Long Range Plan 


Part of Recommendation IV of the 2007 Long Range Plan was the appreciation that “achieving a 
quantitative understanding of the properties of the quark-gluon plasma also requires new investments 
in modeling of heavy-ion collisions, in analytical approaches, and in large-scale computing.” Since then 
there has been tremendous progress along these lines. Convergence has been reached in lattice QCD 
calculations of the temperature for the crossover transition in strongly interacting matter which has now 
been established at 145 MeV < T c < 163 MeV] j33,72,73,2031. Continu um extrapolated results for the 
equation of state, the speed of sound and many other properties of strong interaction matter have also 
been provided 132,331. The modeling of the space-time evolution of heavy-ion collisions has become 
increasingly reliable. (2+l)-dimensional, and subsequently, (3+l)-dimensional relativistic viscous fluid 
dynamics computations have been performed. All such computations use an equation of state extracted 
from lattice QCD. A paradigm shift occurred with the broad appreciation of the importance of fluctuations. 
For heavy-ion collisions at the highest RHIC and at LHC energies, the Color Glass Condensate (CGC) 
effective theory of QCD provides a framework to compute the properties of the fluctuating initial state. 
Under the aegis of the JET topical collaboration, a successful effort was undertaken to consolidate the 
results of different approaches to computing the transport properties of jets as they traverse the strongly 
correlated quark-gluon plasma. Much progress has been made towards a systematic understanding from 
first principles of the properties of strongly interacting matter at non-zero baryon number density. Such 
studies rely heavily on the development of theoretical concepts on critical behavior signaled by conserved 
charge fluctuation 174,88,1111. They are accessible to lattice QCD calculations which opens up the 
possibility, via dynamical modeling, for a systematic comparison of experimental fluctuation observables 
with calculations performed in QCD |115«119|. This will greatly profit from the steady development of 
computational facilities which are soon expected to deliver sustained petaflop/s performance for lattice 
QCD calculations. 

An outcome of these efforts has been the development of a standard dynamical framework for heavy-ion 
collisions, as described Sec. |3.1[ Relativistic viscous fluid dynamics forms the core of this dynamical 


framework and, at high RHIC and LHC energies, it describes the largest part of the evolution history of the 
Little Bang - the explosive expansion of the hot and dense QCD medium formed in the collisions. Much 
recent theoretical work has led to increasingly complete formulations of this theory and improved our 
understanding of its limits of applicability to heavy-ion collisions, and further improvements to optimize 
the framework for such applications are ongoing. 

As described earlier in Sec. |3.1| this standard framework does a good job of describing the wealth of data 
obtained on bulk spectra and event-by-event distributions of anisotropy coefficients, both at RHIC and the 
LHC. To turn this framework into a Little Bang Standard Model requires fixing all the model parameters 
that cannot (yet) be directly computed theoretically, and further refining our descriptions of the initial 
energy deposition and thermalization stages and the interfaces between the different collision stages (see 
following subsection). This process has started but is by no means complete. Finalizing the Little Bang 
Standard Model also requires new sophisticated model/data comparison tools and technologies that are 
just now being developed. Still, with the limited data/theory comparison tools that have so far been 
brought to bear on the large sets of experimental data collected at RHIC and LHC, the specific shear 
viscosity q/s of QCD matter created at RHIC could be constrained to be approximately 50% larger than 
the conjectured lower bound of 1/47T = 0.08, and to be about 2.5 times larger than this bound at the 
LHC. 

Concurrent efforts on extracting the jet quenching parameter q/T 3 from similar theory-data comparisons 
have narrowed the range of values for this parameter to 2 < q/T 3 < 6 within the temperature range probed 
by RHIC and the LHC, nearly an order of magnitude lower than some previous estimates for this quantity. 
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3.5.2 Open questions and future goals 


The development of a broad consensus on the elements of a Standard Model paradigm provides the 
basis for a deeper exploration of each of these elements, as well as the opportunity to further solidify 
the overall picture. Some of the challenging issues over which we need to get better control include i) 
the pre-equilibrium “Glasma” dynamics of coherent gluon fields, and the approach to thermalization; 
ii) the extraction of the values and temperature dependences of transport parameters that reflect the 
many-body QCD dynamics in deconfined matter; iii) the initial conditions at lower collision energies where 
the Glasma framework breaks down; iv) the proper inclusion of the physics of hydrodynamic fluctuations; 
v) an improved treatment of hadron freeze-out and the transition from hydrodynamics to transport theory, 
in particular the treatment of viscous corrections that can influence the extraction from data of the 
physics during the earlier collision stages. 

Quantitative improvements in these aspects of the dynamical modeling of a heavy-ion collision will lead to 
increased precision in the extraction of the underlying many-body QCD physics that governs the various 
collision stages. Additional conceptual advances in our understanding of QCD in matter at extreme 
temperatures and densities are required to answer a number of further outstanding questions. We here 
list a few of them, in chronological order as seen by an observer inside a heavy-ion collision: 

1. A complete quantitative understanding of the properties of the nuclear wave functions that are 
resolved in nucleus-nucleus and proton-nucleus collisions so far remains elusive. Progress requires the 
extension of computations of the energy evolution of these wave functions in the Color Glass Condensate 
(CGC) framework to next-to-leading logarithmic accuracy, and matching these to next-to-leading order 
perturbative QCD computations at large momenta. Simultaneously, conceptual questions regarding the 
factorization and universality of distributions need to be addressed for quantitative progress. These ideas 
will be tested in upcoming proton-nucleus collisions at RHIC and the LHC, and with high precision at a 
future EIC. 

2. How the Glasma thermalizes to the quark-gluon plasma is not well understood. There has been 
significant progress in employing classical statistical methods and kinetic theory to the early stage dynamics 
- however, these rely on extrapolations of weak coupling dynamics to realistically strong couplings. 
Significant insight is also provided from extrapolations in the other direction - to large couplings - using 
the holographic AdS/CFT correspondence between strongly coupled AA = 4 supersymmetric Yang-Mills 
theory in four dimensions and weakly coupled gravity in an AdSsxSs space. Significant numerical and 
analytical progress can be anticipated in this fast evolving field of non-equilibrium non-Abelian plasmas. 

3. A novel development in recent years has been the theoretical study of the possible role of quantum 
anomalies in heavy-ion collisions. A particular example is the Chiral Magnetic Effect (CME), which 
explores the phenomenological consequences of topological transitions in the large magnetic fields created 
at early times in heavy-ion collisions. How the sphaleron transitions that generate topological charge 
occur out of equilibrium is an outstanding question that can be addressed by both weak coupling and 
holographic methods. Further, the effects of these charges can be propagated to late times via anomalous 
hydrodynamics. While there have been hints of the CME in experiments, conventional explanations 
of these exist as well. For a future beam energy scan, quantifying the predictions regarding signatures 
of quantum anomalies is crucial. This requires inclusion of the aromalies into the standard dynamical 
framework. We note that the study of the CME has strong cross-disciplinary appeal, with applications in 
a number of strongly correlated condensed matter systems. 

4. Noteworthy progress has been made in thermal field theory computations of photon and di-lepton 
production in heavy-ion collisions, where NLO computations are now available. The melting of the p 
meson observed in dilepton spectra has been shown to be compatible with chiral symmetry restoration. A 
challenging problem remains the determination of the axial-vector spectral function in the medium, to 
identify the mechanism underlying chiral restoration in the hadron spectrum. The discrepancy between 
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theory and experiment in the “photon V 2 puzzle” mentioned in Sec. 3.1 may require additional sources of 
electromagnetic emissivity, likely in the near-T c region, or arise from more unconventional mechanisms in 
the early collision stages, e.g. transient Bose-Einstein condensates or effects arising from the coupling of 
the conformal anomaly to external magnetic fields. 


5. As observed in Sec. 3.3.1 progress has been made in quantifying the jet quenching parameter q, which 
characterizes an important feature of the transverse response of the quark-gluon medium. Significant 
advances were made in calculating next-to-leading order corrections to q, its QCD evolution and proving 
its factorization from the hard jet production processes [204 -2061. We are beginning to understand the 
flow of the energy radiated by hard partons into a hot medium and its effect on jet shapes 1170,207-2101. 
Soft Collinear Effective Theory (SCET), developed a decade ago in high energy physics, has led to 
improved accuracy of the theoretical predictions of jet cross sections and jet substructure observables 12111. 
It has been generalized to describe jet propagation in nuclear matter 1212,2131 and applied to parton 
transverse momentum broadening 12141. However, significant challenges persist. A general proof of 
factorization of the transport coefficients from the hard partons in an arbitrary medium is still outstanding. 
A second important transport parameter e, characterizing the longitudinal drag of the medium on the 
hard or heavy probe, also needs to be quantified. For confrontation with the available jet fragmentation 
data, the mentioned recent insights into how parton showers develop in the quark-gluon plasma must be 
implemented into Monte-Carlo shower evolution codes that are coupled to a realistic dynamically evolving 
medium. Finally, novel directions for extracting information on the non-perturbative dynamics of the 
strongly correlated quark-gluon plasma from first principles have been pointed out with recent attempts 
to compute the jet quenching parameter using lattice techniques 1215,2161. 

6. Quarkonia and heavy flavor, like jets, are hard probes that provide essential information on the quark- 
gluon plasma on varied length scales. Further, the two probes find common ground in studies of 6-tagged 
and c-tagged jets 12171. Open heavy flavor has attracted considerable theoretical attention 1218-2231, 
but the mechanisms of D- and 73-meson production need to be further elucidated using input from the 
new experimental capabilities at RHIC and LHC. Non-relativistic QCD (NRQCD) has been the standard 
framework for describing heavy quarkonium production in p+p collisions, and much progress has recently 
been made to extend it to p+A 1224, 2251 and A+A collisions 1226, 2271. However, this approach has 
difficulties in describing the polarization of quarkonium produced at collider energies. Recently, a new 
approach based on perturbative QCD factorization 1228-2311 and soft-collinear effective theory 1232,2331, 
called “fragmentation-function approach' ' [183), was developed to address these and other issues. 

Combinations of lattice QCD and effective field theory approaches (NRQCD, pNRQCD and HQET) have 
furthered our understanding of in-medium quarkonium properties [234 -2421 and heavy-flavor diffusion 
coefficients |243,2441. The expected sequential melting pattern has been confirmed 12451, but an accurate 
determination of the in-medium meson properties (masses and widths) will require continued efforts. 
Implementation of quarkonium 1187, 246,2471 and open heavy-flavor transport 1248,2491 into realistic 
bulk evolution models for nucleus-nucleus collisions is ongoing. Tightening the constraints on both the 
theoretical input and the evolution models will require extensive dynamical modeling. 

7. An outstanding intellectual challenge in the field is to map out the QCD phase diagram. While the 
lattice offers an ab initio approach, its successful implementation is beset by the well known sign problem, 
which is also experienced in other branches of physics. Nonetheless, approaches employing reweighting 
and Taylor expansion techniques have become more advanced and are now able to explore the equation of 
state and freeze-out conditions at baryon chemical potentials hb/T < 2. This covers a large part of the 
energy range currently explored in the beam energy scan and suggests that a possible critical endpoint 
may only be found at beam energies less than 20 GeV. Other promising approaches include the complex 
Langevin approach ]250, 251 [ and the integration over a Lefschetz thimble 1252,2531. There has been 
considerable work outlining the phenomenological consequences of a critical point in the phase diagram. 
However, quantitative modeling of how critical fluctuations affect the measured values of the relevant 
observables will require the concerted theoretical effort sketched in Sec. |3.2 
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We cannot emphasize strongly enough that the impressive intellectual achievements outlined and the 
challenges ahead depend strongly on further development of the theory of strongly interacting matter 
which involves advances in heavy ion phenomenology, perturbative QCD, lattice QCD and effective field 
theories for QCD as well as the strong synergy with overlapping and related areas in particle physics, 
condensed matter physics, cold atom physics, string theory and studies of complex dynamical systems. In 
the case of string theory and condensed matter physics, a strong argument can be made that developments 
in heavy-ion collision theory have influenced developments in those fields. The depth and extent of 
interaction of this sub-field of nuclear theory with other branches of physics is perhaps unprecedented in 
nuclear physics. 

3.5.3 What the field needs 

Looking ahead, sustaining and expanding the health of this sub-field of nuclear theory will depend on the 
following key items: 

1. Strong continued support of the core nuclear theory program supporting university Pi’s, national lab 
groups and the National Institute for Nuclear Theory (INT). Together, they play an essential role in 
generating and implementing key ideas driving the field, and in training the next generation of students 
and post-doctoral fellows. 

2. Strong continued support of the DOE Early Career Award (ECA) program in Nuclear Theory, as well as 
the NSF Early Career Development (CAREER) and Presidential Early Career (PECASE) award programs. 
These provide an important opportunity for the field to recognize and promote the careers of the most 
outstanding young nuclear theorists. 

3. Strong support of expanded computational efforts in nuclear physics, as outlined in the Computational 
Nuclear Physics white paper and reflected in Recommendation #3 of this document (see Section [2]). 
Progress in heavy-ion theory is strongly linked to the availability of a diverse and expanding array of 
computational resources, including both leadership class and capacity class computational resources. 

4. Continuation and expansion of Topical Research Collaboration program. These collaborations are 
especially valuable where there are several strains of theory developments that need to be coordinated 
and harnessed to address specific goals. An example of such a successful effort is the JET Topical 
Collaboration involving the co-ordinated effort of both theorists and experimentalists, as discussed earlier 
in this document. The field has several outstanding challenges that require the synthesis of a broad range 
of expertise, and it could strongly benefit from an expansion of such collaborations. 

5. Strong support for the conclusions of the Education and Innovation Town Hall meeting. The numbers 
of scientists from underrepresented minorities and women are particularly low in nuclear theory. We 
encourage efforts to remedy this situation to a level that at least reflects the diversity of talent seen in 
other fields of science. 
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4 Understanding the glue that binds us all: The Electron Ion Collider 

4.1 The next QCD frontier 

Atomic nuclei are built from protons and neutrons, which themselves are composed of quarks that are 
bound together by gluons. Quantum Chromo-Dynamics (QCD), the gauge theory of the strong interaction, 
not only determines the structure of hadrons but also provides the fundamental framework to understand 
the properties and structure of atomic nuclei at all energy scales in the universe. QCD is based on the 
exchange of massless gauge bosons called gluons between the constituents of hadrons, quarks. Without 
gluons there would be no protons, no neutrons, and no atomic nuclei. Matter as we know it would 
not exist. Understanding the interior structure and interactions of nucleons and nuclei in terms of the 
properties and dynamics of the quarks and gluons as dictated by QCD is thus a fundamental and central 
goal of modern nuclear physics. 

Gluons do not carry an electric charge and are thus not directly visible to electrons, photons, and other 
common probes of the structure of matter. An understanding of their role in forming the visible matter in 
the universe thus remains elusive. The Electron Ion Collider (EIC) with its unique capability to collide 
polarized electrons with polarized protons and light ions at unprecedented luminosity, and with heavy 
nuclei at high energy, will be the first precision microscope able to explore how gluons bind quarks to 
form protons and heavier atomic nuclei. 

By precisely imaging gluons and sea quarks inside the proton and nuclei, the EIC will address some of the 
deepest fundamental and puzzling questions nuclear physicists ask: 

• Where are the gluons and sea quarks, and how are their spins, distributed in space and momentum 
inside the nucleon? What is the role of the orbital motion of sea quarks and gluons in building the 
nucleon spin? 

• What happens to the gluon density in nuclei at high energy? Does it saturate? Does this mechanism 
give rise to a universal component of matter in all nuclei, even the proton, when viewed at close to 
the speed of light? 

• How does the nuclear environment affect the distributions of quarks and gluons and their interactions 
in nuclei? How does nuclear matter respond to a fast moving color charge passing through it? How 
do quarks dress themselves to become hadrons? 

A full understanding of QCD, in a regime relevant to the structure and properties of hadrons and nuclei, 
demands a new era at the EIC of precision measurements that can probe them in their full complexity. 
Theoretical advances over the past decade have resulted in the development of a powerful formalism that 
provides quantitative links between such measurements and the above questions physicists are trying 
to answer. A second important advance in recent years is the increasing precision and reach of ab 
initio calculations performed with lattice QCD techniques. Using the experimental data from an EIC, 
physicists will, for the first time, be able to undertake the detailed comparative study between experimental 
measurements and the predictions made by continuum- and lattice-QCD theory, as well as elucidate 
aspects of the structure of hadrons and nuclei whose investigation still requires more phenomenological 
theoretical methods. 

The experimental study of how hadrons and nuclei emerge from the laws of QCD is a global scientific 
priority. Two world-leading facilities in the U.S., CEBAF at Jefferson Lab and RHIC at BNL, are 
international centers for the study of QCD. With the increase of its beam energy to 12 GeV, Jefferson Lab 
will soon operate a unique electron microscope, which will systematically map the structure of protons 
and other nuclei in the valence quark region. In addition to its discovery and continuing exploration of the 
strongly coupled quark gluon plasma (QGP), RHIC has used its unique capability as a polarized proton 
collider to make a first direct determination of the contribution of the gluons to the proton's spin. 
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The high energy, high luminosity polarized EIC will unite and extend the current scientific programs 
at CEBAF and RHIC in dramatic and fundamentally important ways. It will enable us to image the 
transverse momentum and position distributions of quarks and gluons inside fast moving hadrons. When 
hadrons move at nearly the speed of light, the low-momentum gluons contained in their wave functions 
become experimentally accessible. The EIC will study the way in which gluons interact with each other 
by splitting and fusing in addition to providing three-dimensional images of the confined motion of quarks 
and gluons and their spatial distribution. By colliding electrons with heavy nuclei moving with the light 
speed, the EIC will provide access to a so far unconfirmed regime of matter where abundant gluons 
dominate its behavior. Such universal cold gluon matter is an emergent phenomenon of QCD dynamics 
and of high scientific interest and curiosity. Furthermore, its properties and its underlying QCD dynamics 
are critically important for understanding the dynamical origin of the creation of the QGP from colliding 
two relativistic heavy ions, and the QGPs almost perfect liquid behavior. Enabling these research activities 
all in one place, the EIC will be a true “QCD Laboratory”, unique of its kind in the world. 

4.2 Science highlights and deliverables at the EIC 


The high energy, high luminosity polarized EIC will unite and extend the scientific programs at CEBAF, 
RHIC and the LHC in dramatic and fundamentally important ways. 



Figure 10: Left panel: The increase in the proton momentum fraction x vs. the square of the momentum 
transferred by the electron to proton, Q 2 , accessible to the EIC in e+p collisions. Right panel: The 
projected reduction in the uncertainties of the gluons' (AG) and quarks' (AS) contributions to the 
proton's spin. 


Proton Spin: Recent measurements at RHIC along with state-of-the-art perturbative QCD analyses 
have shown that gluons carry approximately 20-30% of the proton’s helicity, similar to the contribution 
from quarks and anti-quarks. The blue band in the right panel of Figure K) shows the current level of 
uncertainties. The knowledge is limited by the x-range explored so far. The EIC would greatly increase 
the kinematic coverage in x and Q 2 , as shown in the left panel of Figure 10 and hence reduce this 
uncertainty dramatically, to the level depicted by the red and yellow bands in the right panel. 


Motion of quarks and gluons in a proton: Semi-inclusive measurements with polarized proton beams 
would enable us to selectively and precisely probe the correlations between the spin of a fast moving 
proton and the confined transverse motion of both the quarks and gluons within it. Images in momentum 
space as shown in the left panel of Figure |TT] are simply unattainable without the polarized electron and 
proton beams of the proposed EIC. 
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Figure 11: Left panel: Transverse momentum distribution of a sea tt-quark with longitudinal momentum 
fraction x = 0.1 in a transversely polarized proton moving in the z-direction, while polarized in ^-direction. 
The color code indicates the probability of finding the up quarks, with dark red indicating the highest 
probability. Right panel: Projected precision of the transverse spatial distribution of gluons obtained 
from exclusive J/ij) production at the EIC. 



Figure 12: Left panel: Schematic landscape of probe resolution vs. energy, indicating regions of non- 
perturbative and perturbative QCD, including low to high parton density and the transition region. Right 
panel: Ratio of diffractive over total cross section for deep-inelastic scattering (DIS) of electrons on gold, 
normalized to DIS on the proton, for different values of the square of the invariant mass of the hadrons 
produced in the collisions, with and without saturation. 


Tomographic images of the proton: By choosing particular final states in e+p scattering, the EIC, 
with its unprecedented luminosity and detector coverage, will create detailed images of the proton’s gluon 
matter distribution, as shown in the right panel of Figure [Tl] Such measurements would reveal aspects of 
proton structure that are intimately connected with QCD dynamics at large distances. 


QCD matter at extreme gluon density: When fast moving hadrons are probed at high energy, the 
low-momentum gluons contained in their wave functions become experimentally accessible. By colliding 
electrons with heavy nuclei moving at light-speed, the EIC will provide access to a so far unconfirmed 


regime of matter where abundant gluons dominate its behavior as shown in the left panel of Figure 12 


Such cold gluon matter is an emergent phenomenon of QCD dynamics and of high scientific interest and 
curiosity. Furthermore, its underlying QCD dynamics and its predicted universal properties are critically 


34 



































important for understanding the dynamical origin of the creation of the QGP from colliding two relativistic 
heavy ions. 

By measuring diffractive cross-sections together with the total deep-inelastic scattering (DIS) cross- 
sections in e+p and e+A collisions, shown in the right panel of Figure [T2j the EIC would provide the first 
unambiguous evidence for this novel state of saturated gluon matter in QCD. The planned EIC is capable 
of exploring with precision the new field of collective dynamics of saturated gluons at high energies. 


Hadronization and energy loss: The mechanism by which colored partons pass through colored media, 
both cold nuclei and hot QGP matter, and color-singlet hadrons finally emerge from the colored partons 
is not understood. A nucleus in the EIC would provide an invaluable femtometer filter to explore and 
expose how colored partons interact and hadronize in nuclear matter, as illustrated in the left panel of 
Figure 13 By measuring ir and D° meson production in both e+p and e+A collisions, the EIC would 


provide the first measurement of the quark mass dependence of the response of nuclear matter to a fast 
moving quark. The dramatic difference between them, shown in the right panel of Figure 33 would be 
readily discernible. The color bands reflect the limitations on our current knowledge of hadronization - 
the emergence of a pion from a colored quark. Enabling all such studies in one place, the EIC will be a 
true QCD Laboratory, a unique facility in the world. 


4.3 EIC machine parameters and designs 

Two independent designs for a future EIC have evolved over the past few years. Both use existing 
infrastructure and facilities available to the US nuclear scientists. At Brookhaven National Laboratory 
(BNL), the eRFHIC concept adds a new electron beam facility, based on an Energy Recovery Linac (ERL) 
to be built inside the RHIC tunnel, to collide electrons with one of the existing RHIC beams. At Jefferson 
Laboratory the Medium Energy Electron Ion Collider (MEIC) concept envisions a new electron and ion 
collider ring complex, together with the 12 GeV upgraded CEBAF, in order to achieve similar collision 
parameters. The machine designs aim to reach the following goals and parameters: 




h 



Figure 13: Left panel: A schematic illustrating the interaction of a parton moving through cold nuclear 
matter, where the hadron is formed outside (top) or inside the nucleus (bottom). Right panel: The 
ratio of the semi-inclusive cross sections for producing a pion (light quarks, red) and D° mesons (heavy 
quarks, blue) in e+lead collisions to e+deuteron collisions, plotted as a function of the ratio z of the 
momentum carried by the produced hadron to that of a virtual photon. 
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• Polarized (~ 70%) beams of electrons, protons and light nuclei; 

• Ion beams from deuteron to the heaviest nuclei (uranium or lead); 

• Variable center of mass energies from ~ 20 to ~ 100 GeV, upgradable to ~ 140 GeV; 

• High collision luminosity ~ 10 33 — 10 34 cm _2 sec _1 ; 

• Capacity to have more than one interaction region. 

4.4 Why now? 

Today, a set of compelling physics questions related to role of gluons in QCD has been formulated, 
and a corresponding set of measurements at the EIC identified. A powerful formalism that connects 
those measurements to the QCD structure of hadrons and nuclei has been developed. The EIC was 
designated in the 2007 Nuclear Physics Long Range Plan as “embodying the vision for reaching the next 
QCD frontier In 2013 the NSAC Subcommittee report on Future Scientific Facilities declared an EIC 
to be “absolutely essential in its ability to contribute to the world-leading science in the next decade." 
Accelerator technology has recently developed so that an EIC with the versatile range of kinematics, 
beam species and polarization, crucial to addressing the most central questions in QCD, can now be 
constructed at an affordable cost. Realizing the EIC will be essential to maintain U.S. leadership in the 
important fields of nuclear physics and accelerator science. 
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Phases of QCD Matter” meets separately 


6 Appendix: Agenda of the “Phases of QCD Matter” and joint sessions 
at the QCD Town Meeting 


Temple University QCD Town Meeting — Saturday, September 13, 2014 


8:00 AM 

Coffee, Tea, Bagels, Doughnuts and Fruit are Available 

9:00 AM 

Welcome and brief introduction to the format and goal of the meeting 

Ulrich Heinz 

Theory Overview 


9:10 AM 

-- Open questions in theory, why they are interesting, what tools we have to 

Krishna Rajagopal 



Discussion (10 min) 

9:50 AM 

-- Lattice QCD for RHIC and LHC: from now to the next decade 

Swagato Mukherjee 



Discussion (lOmin) 

10:30 AM 

Coffee Break 

QGP evolution and progress on EOS and transport properties 


11:00 AM 

-- The standard model for QGP evolution: theoretical status and future 

Bjoern Schenke 



Discussion (10 min) 

11:40 AM 

-- The standard model for QGP evolution: experimental status and future 

Raimond Snellings 



Discussion (10 min) 

12:20 PM 

Lunch Break 

BES and QCD Phase Diagram 


1:20 PM 

-- BES 1 Results, Motivation for BES II, Future Facilities 

Dan Cebra 



Discussion (10 min) 

2:00 AM 

-- Studying the QCD Phase Diagram via BES 

Misha Stephanov 

2:20 AM 

-- Dynamical modeling at low collision energy and high mu_b 

Hannah Petersen 



Discussion (10 min) 

2:50 PM 

Community Input (hot QCD related) 

3:30 PM Coffee Break 


4:00 - 7:30 PM Joint Session 


Figure 14: Agenda of the “Phases of QCD Matter" parallel session, part I (Saturday) 
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Saturday 


Joint Hadronic and Phases of QCD Session 

4:00 PM 

— Nuclear theory for the next decade: an overview 

David Kaplan 



Discussion (10 min) 

4:40 PM 

— Highlights in cold QCD since 2007 

Rolf Ent 

5:00 PM 

— Cold QCD: the next decade 

Naomi Makins 



Discussion (10 min) 

5:30 PM 

— RHIC and LHC overview (“Where we are, where do we need to go?") 

Bill Zajc 



Discussion (10 min) 

6:20 PM 

— Why we need an EIC: a view from 30,000 feet 

Berndt Mueller 

6:40 PM 

— Why we need an EIC: a view from 10,000 meters 

Bob McKeown 

7:00 PM 

Combined discussion period for the preceding 2 talks 

Discussion (30 min) 

Sunday 


Morning EIC Session 

8:30 AM 

Central issues in nucleon structure (theory) 

Xiangdong Ji 

9:00 AM 

Probing nucleon structure at an EIC 

Zein-Eddine Meziani 

9:30 AM 

Probing the properties of QCD with atomic nuclei: theoretical perspective 

Yuri Kovchegov 

10:00 AM 

Probing the properties of QCD with atomic nuclei: experimental perspective 

Thomas Ullrich 

10:30 AM 

Combined discussion period for the preceding 4 talks Discussion (20 min) 


10:50 AM 

Coffee Break 

11:15 AM 

Community Input (EIC related) Weiss/Hyde/Avakian/Jacak 


12:30 PM 

Lunch Break 

1:30 PM 

Discussion about EIC White Paper Abhay Deshpande 


Reports 

2:45 PM 

Report from the Computational Nuclear Physics Town Meeting 1 

Peter Petreczky 

3:00 PM 

Report from the Computational Nuclear Physics Town Meeting 2 

Martin Savage 

3:15 PM 

Report from Educational Town Meeting 

Thia Keppel 

3:30 PM 

Combined discussion period for the preceding reports 

Discussion (30 min) 

4:00 PM Coffee Break 


Figure 15: 


Agenda of the joint session (Saturday/Sunday) 
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Sunday, September 14, 2014 


8:30 AM - 4:00 PM Joint Session 


4:00 PM 

Coffee Break 


"Phases of QCD" session on low x, p+A, and spin 


4:30 PM 

-- Theory: Unique opportunities in p+A collisions at RHIC and LHC 

Zhong-Bo Kang 




Discussion (10 min) 


5:00 PM 

-- p+A experiments at RHIC and the LHC and related upgrades 

John Lajoie 




Discussion (10 min) 


5:30 PM 

-- Opportunities for RHIC Spin in the Coming Decade 

Elke Aschenauer 




Discussion (10 min) 

Monday, September 15, 2014 


Jet Studies of QCD matter 



8:30 AM 

-- Jets at RHIC & LHC: probing the QGP at multiple length scales (theory) 

Xin-Nian Wang 




Discussion (10 min) 

o 

if) 

if) 

Z3 

9:00 AM 

-- Jet studies at LHC & RHIC now and tomorrow 

Gunther Roland 



Discussion (10 min) 

O 

if) 

9:30 AM 

-- "What jets studies at RHIC can tell us about QGP properties" 

Jamie Nagle 

o 



Discussion (10 min) 

“O 

c 

10:00 AM 

Coffee Break 

03 

if) 

Heavy quarks and quarkonium as a probe of the QGP 
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Figure 16: Agenda of the “Phases of QCD Matter” parallel session, part II (Sunday/Monday) 
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